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Abstract 


The  original  motivation  for  this  study  comes  from  the 
climatological  observation  of  three  semipermanent  sub- 
tropical highs  in  each  of  the  Northern  and  Southern  hemispheres 
and  its  explanation  in  the  work  of  H.  J.  Stewart,  which  was 
extended  by  G.  K.  Morikawa,  E.  Swenson  and  A.  S.  Peters. 

A geos trophic  approximation  has  been  formulated  by  A.  S. 
Peters,  for  the  one  layer,  hydrostatic  atmosphere,  over  a 
rotating  sphere.  In  this  yeostrophic  model,  the  coriolis 
parameter  is  treated  as  a constant;  the  height  of  the  layer 
and  its  tangential  velocity  components  are  expressed  in  terms 
of  a stream  function. 

Singular  solutions  composed  of  a finite  number  of  vortices 
that  are  in  equilibrium,  may  be  found.  The  stability  under 

perturbation  of  the  following  four  types  of  configurations 

/ 

of  "geos trophic"  vortices  isstudied: 

(a)  N vortices  of  equal  strength  y equally  spaced  on  a circle 
of  colatitude  0; 

(b)  N vortices  of  equal  strength  y equally  space  on  a circle 
of  colatitude  0,  one  vortex  with  strength  yQ  at  the  pole; 

(c)  the  same  distribution  as  in  (a)  but  satisfying  the  boundary 
condition:  vanishing  normal  velocity  component  at  the 
equator; 

(d)  the  same  distribution  as  (b)  but  satisfying  the  boundary 
condition:  vanishing  normal  velocity  component  at  the 
equator. 
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In  this  study  we  can  distinguish  two  main  parts: 

(1)  The  solution  of  an  algebraic  eigenvalue;  problem  that 
arises  from  a linearized  stability  analysis  of  the 
nonlinear  kinematic  equations  of  vortex  motion; 

(ii)  The  numerical  solution  of  the  nonlinear  equations 

of  vortex  motion,  which  confirms,  in  some  cases,  the 
validity  of  the  linearized  analysis  and  yields  the 
magnitude  of  the  perturbed  motions. 


0 . Introduction 


The  study  of  geostrophic  vortices  was  begun  by 
H.  J.  Stewart  [1]  in  an  attempt  to  explain  the  occurrence 
of  the  three  so-called  semi -permanent  high  pressure  cells, 
in  each  of  the  Northern  and  Southern  hemispheres.  In  this 
work,  Stewart  represented  the  earth  by  a rotating  tangent 
plane  and  the  large  scale  closed  isobaric  systems  of  the 
atmospheric  layer  by  discrete  rectilinear  vortices,  which 
are  defined  by  the  Bessel  function  ,'q.  Later  on,  G.  K. 
Morikawa  and  E.  V.  Swe..son  (2]  studied  in  greater  detail 
the  stability  of  these  single  layer  rectilinear  geostrophic 
vortices  defined  on  a rotating  plane  with  constant  Coriolis 
parameter.  The  rectilinear  vortices  simulated  the  large 
scale  motions  of  the  atmosphere  over  the  earth  by  means  of 
projection  from  a rotating  tangent  plane. 

A.  S.  Peters  [3]  found,  upon  neglecting  the  variation 
of  the  Coriolis  force  with  colatitude,  that  a geostrophic 
vortex  on  a rotating  sphere  is  given  by  a singular  spherical 
harmonic  of  degree  v and  order  0.  A brief  description  of 
this  derivation  {3]  follows: 

Let  S be  a sphere  of  radius  a,  that  rotates  with 
constant  angular  velocity  oj  about  a polar  axis.  Let  p be 
the  distance  from  any  point  in  space  to  the  center  of  the 
sphere  S and  0,4)  be  the  colatitude  and  longitude  of  a 
point  on  the  surface  of  S. 
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Let  p = a and  p = a + h(<J>,6,t)  represent  the  lower 
and  upper  surfaces  of  an  incompressible,  inviscid  fluid 
layer  that  is  gravitationally  attracted  by  S and  suppose 
that  h(<J>,0,t)  <<  a. 

At  time  t = 0,  the  constant  rotary  motion  of  the 
fluid  is  disturbed  by  the  creation  of  concentrated  vortices 
with  axes  that  are  normal  to  the  surface  of  S. 

The  velocity  components  of  a fluid  particle,  relative 
to  S,  are  defined  by 

u = P sin  0 (longitudinal  component  to  the  East) 

v = -p  (latitudinal  component  to  the  North) 

w = ir-  (radial  component), 

at 

The  basic  Eulerian  hydrodynamical  equations  defining 
the  motion  consist  of  a continuity  equation  and  three 
momentum  equations . 

The  shallow  water,  or  long  wave,  model  is  derived  from 
the  Eulerian  equations  by  assuming  that: 

(a)  The  only  body  force  acting  is  that  due  to  the 
gravitational  potential  G of  S and  furthermore 
G = gp;  g = gravitational  constant; 

(b)  The  radial  momentum  equation  can  be  replaced  by  the 
hydrostatic  law,  for  the  pressure  p,  in  terms  of 
the  density  a per  unit  surface  area,  that  is 

p($,e,p,t)  =0 • g • (h+a-p)  , 
which  implies 


-2- 


p = 0 at  the  free  surface,  p = h+a; 

(c)  the  motion  is  such  that  the  nonlinear  terms  in  the 
momentum  equations  can  be  neglected; 

(d)  the  radial  velocity  is  negligible  and  the  radial 
variation  of  u and  v can  be  ignored. 


Under  these  considerations,  an  approximation  to  the 
motion  of  the  atmospheric  layer  on  the  surface  of  S is 
determined  by  a continuity  equation  and  two  momentum 
equations  [ 3] : 


(0.1) 


nt  a sin  6 


< 

» 


9 ( v sin  6)  _ 9u 

96  " 9<f> 


(0.2) 


9u 
9 1 


- 2 cuv  cos  6 = - 


g 9h 
a sin  6 d<p  ' 


(0.3)  f£  + 2„u  cos  e = a ||  , 

whe  re 

h f <£  , 0 , t)  -hQ 

n = r- , hn  = h (4>,t>,0)  = constant  . 

n0  u 

A.  S.  Peters  (3]  has  shown  that  by  neglecting  the 
variation  of  the  coriolis  parameter,  f = 2 co  cos  0 with 
colatitude,  equations  (0.1),  (0.2),  (0.3)  lead  to 


(0.4) 


sin  6 90 


„ 9n  . 1 9 n a , , 2 . 9 n , 

sin  0 + K ^ - ~r—  [4aj.n  + -k  ) 

90  m2  ghp  1 ^ ^2 


sin  6 9 


ifi 

9t* 


2 W1 


= 2a  Z(<p,6,0)  , 

9“0 


where  <d^  = w *cos  6^  - f/2  = oj  »c  = constant  and 
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does  not  depend  on  tine  and  in  fact  is  the  stream  function, 
so  that  the  associated  steady  velocity  components  are 

1 3U/ 

u " a 36 

1 91 

v a sin  6 9 

Equation  (0.5)  describes  the  geostrophic  motion  of  the 

atmospheric  layer  over  a sphere. 

It  can  be  verified  by  direct  substitution  that  (0.5) 

is  satisfied  by 
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(0.6) 


i p(<p,Q,<p.  ,6.)  = r~ P„  (-cos  a.)  , 

jr  j 4 sin  v tt  v'  1' 

where  v is  complex, 


1 

a,  = — , d.  = geodesic  distance  from  the  center 

la  1 

of  the  vortex  to  any  point  (<f>,0)  on 
the  su> jace  of  S. 

Py  (x)  = spherical  harmonic  of  order  v [3]. 

A superposition  of  vortices,  with  centers  at 
(^k  (t)  ' 0Jc  (t) ) satisfi-es  the  system  of  ordinary  differential 
equations  (1.3)  and  (1.4). 

The  stability  of  four  configurations  of  these  geostrophic 
vortices  on  the  sphere  is  studied  in  this  work: 

(a)  N vortices  of  equal  strength  u equally  spaced  on  a 
circle  of  colatitude  6; 

(b)  N vortices  of  equal  strengh  v equally  spaced  on  a circle 
of  colatitude  0,  one  vortex  with  strength  uQ  at  the  pole; 

(c)  the  same  distribution  as  in  (a)  but  satisfying  the 
boundary  condition:  vanishing  normal  velocity  component 
at  the  equator; 

(d)  the  same  distribution  as  (b)  but  satisfying  the  boundary 
condition:  vanishing  normal  velocity  at  the  equator. 
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For  the  problems  under  consideration  there  are  four 


paramete rs : 

1.  Initial  colatitude  of  the  belt  of  vortices  (0), 

2.  Number  of  vortices  on  belt  (N) 

3.  Coriolis  parameter  c = ~- 

c p 2t*> 

4.  Strengths  y and  y^. 

The  linearized  stability  of  an  equilibrium  configuration 

is  studied  for  N>2,  5°  < 6 < 90°,  .2  < c < 1.0  and 

- - - P - 

“°°  < vq/v  < °° * Zones  of  linear  stability  in  the  (cp,6)  plane 
and  the  (yQ/y,  6)  plane  are  given  in  the  Appendix. 

Numerical  integration  of  the  nonlinear  equations  of 
vortex  motion  (1.3)  and  (1.4)  is  used  to  test  the  validity 
of  the  linear  analysis,  i.e.,  to  see  when  nonlinear  effects 
are  important,  and  to  determine  the  linear  and  nonlinear 
motions.  Results  on  the  nonlinear  initial  value  problem 
are  given  in  Section  3 and  in  the  Appendix. 

These  computations  show  that  nonlinear  effects  are 
important  near  the  boundary  of  the  linear  stability  zone 
and  that  the  nonlinear  effects  are  more  pronounced  when  a 
polar  vortex  is  considered. 
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1.  Rations  of  Motion  for  Four  Configurations  of 

In  general,  the  position  of  each  vortex  on  the  surface 
of  S is  described  by  its  colatitude  6.  and  longitude  * 
The  stream  function  iM4>,6)  is  related  to  the  velocity 
components  u,v  (longitudinal,  latitudinal)  by 


(l.i)  ( u , v)  = (I  if  , - iil  - ( 

a 30  a sin  0 3ij|J  ~ ( 


!a  sin  6<j>,  -ad) 


The  stream  function  generated  by  the  superposition  of  N+l 
vortices  must  satisfy 


2 4a  iiij 


l L e,'n  ^ , 1 d 

sin  A 30  Sln  3 90  + ~~r~  2 SiT 

sin  0 3 <p  gn0 


N 6 (0-4>  ) 6(6-0.) 

= > Ll  . — — 

i sin  0 . 

l 


i=0 


As  we  can  see  from  (0.6),  the  above  equation  is  satisfied  by 


(1.2) 


i M4>,0)  = 


4a  sin 


V Tf 


N 

l 

i=0 


V : PX)  (-cos  a.  ) 


With  cos  ou  = cos  0 cos  + sin  0 sin  0^^  cos  ($-$  ) , 
where  (<p,d)  is  an  arbitrary  point  on  the  sufrace  of  S 
and  ( <fo  , 6^ ) are  the  coordinates  of  the  center  of  the  ith 

vortex,  0 < i < N;  the  constant  is  called  the  strength 
of  the  ith  vortex. 


It  can  be  shown  [3]  that  the  nonlinear  kinematic  equations 
of  motion  of  the  center  of  the  kth  vortex  can  be  written,  by 
using  (1.1)  and  (1.2),  as  if  the  vortex  particle  (<}>  ,0  ) moves 

K K 
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under  the  influence  of  the  remaining  vortices: 


(1.3)  a sin  0kik  = v.,  y^-  — Pi  V'cos  «ik) 

•(sin  0k  cos  6^-sin  0^  cos  ©k  cos  (<J>k-<fK  ) ) , 
. N 

(1‘4)  "a\  ■ 4a  sin  77  J0  pi  pv('cos  “ik>  sin  Vsin 

i /*. 

with 

(1.5)  cos  a-k  = cos  9k  cos  0i  + sin  0k  sin  0.^  cos  (4>k“<l>i) 


for  k = 0,1,..., N,  and 


Pu(x) 


Sp«“ 


1.1.  Cases  (a.b) . N Vortices  on  a Circle  of  Colatitude, 
with  and  without  a Polar  Vortex. 

Consider  the  configuration  that  consists  initially  of 
one  polar  vortex  of  strength  Ug  and  N vortices  of  equal 
strength  y,  equally  spaced  on  a circle  of  colatitude  y, 
that  is  let 

6k  = Y , 4>k  = (k-1)  — jj-  , V*k  = ^ = 1,  • • • »N) 

(1.6) 

0g  = 0 , fp0  = arbitrary,  yQ  = arbitrary. 

Inserting  (1.6)  into  (1.3),  (1.4),  (1.5)  we  find 
that  the  latitudinal  component  of  the  velocity  vQ  at 
the  pole  is 
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sin  = 0 , 


N 


V0  4a  sin  vtt 
since 


l p' (-cos  y)  sin  Y 


i=  1 


N 


l sin  (4-n“4>i-)  = 0 
k-1  0 k 


N 


itudiial  component  un  also  vanishes  since  ^cos  ^i' 


The  longi luui  un  uuuujuhcih. 

For  the  circle  vortices,  we  find  that 

N 


v,.  = 


k 4a  sin  vtt 


by  symmetry,  and 


l (pki)  sin  y sin  wki  = 0 


k=l 


u.  = a sin  Y 4k  = 


4a  s 


i |unP'  (-cos  y) 

in  vtt  ^ 0 v 


sin  y 


N 


where 


+ u l p' (Pki)  sin  Y cos  Y (1-cos  wki) 
i=  1 
i^k 


. 2 


ki 


= -cos  Y " sin  Y cos  wki 


,,  . . 2 TT 

wki  k-i  N 

Thus,  the  vortices  remain  on  the  circle  of  colatitude 
Y but  they  move  around  the  circle  with  constant  angular 
velocity 

1 


(1.7)  4>k  — i 2 


4a  sin  vtt 


Mq  Pv(-cos  y) 


N l 

+ M pv(pki}  cos  Y (1_cos  wki)J 


i/k 
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1 • 2 • Cases  (c,d)  . N Vortices  on  a Circle  of  Colatitude, 

with  and  without  a Polar  Vortex;  Satisfying  a Boundary 
Condition  at  the  Equator. 

By  means  of  reflection  across  the  equator  we  can  look  for 
solutions  that  satisfy  some  prescribed  boundary  condition  at 
0 = 90°. 

The  boundary  condition  treated  here  is 


2k 


| 


0=tt/2 


0 


i.e.  vanishing  latitudinal  velocity  component  at 
the  equator,  h = tt/2.  The  stream  function  for  the  case 
of  a polar  vortex  along  with  N other  vortices  in  the  northern 
hemisphere  satisfying  the  boundary  condition  is  given  by 


(1.&) 


cos 

■ 

cos  Y • 

l 


1 * 

4a  sin  vrr  JQ  MPv(-COS  Yi)-Pv(cos  Y • ) ) 

cos  0 cos  9.  + sin  0 sin  0.  cos  (*-<!>.) 

i l l 

cos  0 cos  Ch  - sin  0 sin  9^  cos  ( ^ ) . 


This  sum  can  be  interpreted  as  representing  a superposition 
of  N+ 1 symmetrically  placed  vortices  in  the  southern  hemi- 
sphere with  strengths  that  are  the  negatives  of  those  of 
the  corresponding  vortices  in  the  northern  hemisphere. 

Equations  (1.1)  and  (1.8)  yield  the  nonlinear  equations 
of  motion  of  the  center  of  the  kth  vortex, 
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COS  0 . 

1 


(1.9)  a sin  9^0^  = 


4a  sin  vtt  { ^^i  [(Pv  ^ -cos  aik^+Pv^cos  a^))sin  0j 


i/k 


(pv(-cos  aik)-Pv(cos  aik))cos  0k  sin  0..  cos  (4>k-0i )J 


+ ^Mk  Pv(cos  20k)  sin  0k  cos  ©k| 

1’1°)  3 k 4a  sin  vtt  ^ jpv(-cos  aik)_Pv(cos  aik)j 


i^k 


•sin  0 . sin  (<}>,  -<fc . ) 
i k Yi' 


(I’lD  COS  aik  ~ cos  ek  cos  9i  " sin  0k  sin  0i  cos  (4-k-4-i) 


For  the  case  that  the  N nonpolar  vortices  have  equal  strengths 

U and  lnitially  are  on  the  line  of  colatitude  y,  (1.6)  again 
yields  = 0 , 1 £ k < N,  but  now 


(1.12)  *k  = « = —1 


4a  sin  vtt  l 


.{uo(pv<- 


cos  Y)+P  (cos  y) 


u 


+ U cos  y 


/ N . n . 

iJi  Pv(pki)(1-cos  Wki)+J1Pv(pki)(1+cos  wk.)] 


'1=1 

k=l 

i^k 

2 

2 

’ki  = 

-cos  y - 

sm  y cos  w,  . , 
ki  * 

* 

2 

2 

ki  = 

COS  Y - 

sm  y cos  w,  . , 
ki  ' 

n 

•H 

Oc-i,  21 

■ 

That  is,  the  vortices  remain  on  the  circle  of  colatitude 
Y and  they  traverse  it  with  the  constant  angular  velocity  (2 


given  by  (1.12) 
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1 . 3 Invari an ts 


The  kinematic  equations  of  motion  can  be  written  in 
a more  systematic  way,  by  using  the  Kirchhoff  function  W 
for  geostroohic  vortices  on  a sphere. 

For  a configuration  of  (Nfl)  vortices  and  no  boundary 
condition  at  0 = tt/2  , the  Kirchhoff  function  is  defined  by 


i N 

( 1 • 13)  W = — : 7 ij  ’ > p (-COS  fl 

8a  sm  v 7T  . 4 . MiH-i  v'  cos 

i/3=0  - v 1 


cos  0 


D 


-sin  ei  sin  0.  cos  (4^-$.)) 


Then,  with  uR  - a sin  0^  and  vR  = -aOk  , it  follows  that 


^ _ riw  1 3W 

*(W  ~ ' sIrTl£  3^)  ' sin  0k  ^ 0 . 


Since  W depends  only  on  the  relative  distances  between 
the  vortices  and  is  not  an  explicit  function  of  time,  it 
is  easy  to  verify  that 


dt  “ 0 and  therefore  W is  itself  an  invariant,  i.e. 
W = constant; 

(ii)  and  from 


■ "k  a 8k 


we  have 


1 3W 

nn;  \ * ° 


a IB  Jk0k  sin  8k  « - f |W  0 , 


k=0 


k=0 


as  an  easy  consequence  of  the  way  that  the  angles  $ 
appear  in  W in  (1.13).  Therefore  by  integrating,  we  get 
another  invariant  defined  by 
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constant. 


N 

d.i4)  f1(e1,e2, . . . ,ek)  = l Mk  cos  ek  = 

k — 0 

Clearly  f^  is  an  invariant  of  the  motion, 
(iii)  Multiplying  both  sides  of 


= a sin  6k^  by 


l-cos  ek 
sin  0k  ' 


sin  8k  / 0 , and 


summing  over  k we  define 

(i.i5)  f2(e1#...#ek, 


'V  ^ 


? 9W_ 

1 dd, 


k=0 


l-COS  0^ 

sin  0, 


N 


= I 

k=0 


»k(1 


cos  0k) 


For  the  case  with  boundary  condition  we  can  derive 
analogous  expressions  fcr  the  functions  W,  f^  and  f 2 . 

We  shall  see  that  although  the  function 
f2  (0^, . . . , 0k,  <£lf...,<£k)  given  by  (1.15)  is  not  an 
invariant,  its  linearization  with  respect  to  the  equili- 
brium solution  (2.1),  (2.2),  is  an  invariant  for  the 
linearized  motion. 

The  equations  of  motion  state  that  the  velocity  of 
the  center  of  any  vortex,  say  the  kth,  is  the  sum  of  the 
field  velocities  of  the  remaining  ones  evaluated  at  the 
position  of  the  kth  vortex. 

If  the  strength  of  a vortex  is  positive  the  velocity 
field  of  the  vortex  has  counterclockwise  or  cyclonic 
rotation . 

In  this  system,  each  vortex  center  moves  like  a massless 
particle  in  the  velocity  field  of  the  others. 
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2 • Linear  Stability  Analysis  of  the  Nonlinear  Equations  of 


Motion 


The  linearized  equations  are  obtained  by  makinq  a 
perturbation  expansion  of  (t^/O^)  about  the  equilibrium 
solution  satisfying  the  initial  condition  (1.6). 

That  is,  for  the  kth  vortex 


(2.1) 


°k  = y * c3k  ■ 


” <k-l)  ^ + Sit  + Cdk  ; 


while  for  the  polar  vortex 


H 0 c B 0 ' 


<t>Q  = ca0  + Qt  ; 


where  c is  small  and  and  Bk  are  functions  of  the  time  t, 


2.1.  Cases  (a,b) . N+l  Vortices,  No  Boundary  Condition 


at  Equator. 


The  formulas  (2.1),  (2.2)  are  inserted  into  (1.3),  (1.4), 

(1.5).  By  expanding  in  powers  of  e and  retaining  first 
order  terms,  we  get  the  following  linearized  equations 
of  motion: 
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2 ' / 

(2.3)  -4a  sin  vn6k  = uQ  Pv(-cos  y)H  sin  wkl~n  cos 


N ft  N ft 

+ .1  <i8i  + I Bu,ak-ai» 
1=1  1=1 


1 < k < N , 


wkl} 


i^k 


i/k 


2 . 


(2.4)  4a  sin  vn  sin 


YV‘l,o{%(_ 


cos  y)cos  y 


■}{« 


+ ( -cos  y)  sin  YfjS  cos  sin 


wki} 


N 


N 


+ Akk6k  + J,  ^i  - J,  Rkiai  ' 1 i k ^ N ' 

i^k  i/k 


(2.5) 


4a  sin  vtt  . • 


( n + HO  * (Pv(-cos  y)  N cos  y 


+ P (-cos  y)  sin' 


> r 2 s\ 

Y l cos  wilJ|  • 5 

, f N N 

+ P ^ ( —cos  y)  j^sin  y I ci  sin  w^-  cos  y Bicos  wi;L 


. 2 r 

- Pv(-cos  y)  sin  Y z Bi  cos  Wil  ' 

i=l 


N 


(2.6) 


4a  sin  vtt 


. . » 

K-an)  = -(py(-cos  y)  n cos  y 


N 


+ P (-cos  y)  sin^Y  I sin  w^i  ) *n 


.2  ” _ . _ 2 

‘ Y 

i=l 
N 


N 


+ P,(-cos  y)  (sin  Y [ a.  cos  w. . + cos  Y Z B i sin  w..) 
v i=1  i i-L  i=! 


2 _ 

+ P (-cos  y)  sin  Y I cos  w.  . ; 

v i=1 


N 
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where 


j i 

(2.7)  Aki  = p p ( P^i)  cos  y sin  wki 


+ p P (p,  .)  sin  y cos  Y sin  w.  . (1-cos  w,  . ) , 
v ki  ' ki  ki 


Bki  * “ pj(pki)  sin  Y cos  “ki"  uPv(pki)  sir,3Y  sin2"ki 


* 2 
A A = p„(cos  y P (-cos  Y ) + sin  y P (-cos  y)1 


'kk  MC 
N 


i 2 2 

+ p l P (Pk^)  (cos  y + sin  Y cos  wk . ) 
i = l 
i^k 

N 


-4azsin  vtt  cos  ytt  + 
ki' 


”22  2 
+ g l P (Pki)  sin  Y cos  y (1-cos  wki)  , 


(2.8)  Bk. 


= p P (p,  .)  cos  y sin  y sin  w,  . 
v ki  ki 


"3  . 

+ p IJ  (pk^)  sin  y cos  Y (1-cos  wk^)  sin  wk^  , 


• 2 2 
Aki  = -uP  (Pki)  (sin  y + cos  Y cos  wRi) 


"22  2 
+ g P (pki)  sin  Y cos  Y (1-cos  wki)  , 


and 


£ = 0Q  cos  aQ  , n = 80  sin  aQ  . 
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2.2.  i d s e s — (c,d)  . — N+l  Vortices,  Satisfying  Boundary  Condi- 
tion at  the  Equator. 


By  inserting  (2.1)  ard  (2.2)  into  equations  (1.9), 
(1.10),  (1.11)  and  linearizing  analogously,  it  follows  that 

(2.9)  -4a2  sin  vtt  ek  = p0|p^(-cos  y)V(cos  y)J 

}N  N 

' lik-N  ' 


sin  wkl-n  cos 


i^k 


2 C i 

(2.10)  4a  sin  vtt  sin  = -P0|P^(-cos  y)V(cos  y)J 


•cos  yjr  cos  wkl+n  sin  wk]J 


- v 


{Py ( ~ cos  Y)+Py(cos  y)|sin2y  j^cos  «kl+n  sin  «kl| 
+ AkkBk  + il1  AkiBi  ‘ Bkiai  ' 1 i k < N , 


i^k 


i^k 


2 

n m 4a  sin  vtt  . • _ _ 
(2.H)  7, (n+ns) 


- {(v- 


cos  y ) 


+ Pv(cos  Y) )n  cos  y + (p^ (“cos  y)-P^(cos  y)) 


. 2 
sin 


Y cos2  wu  + 2 r v1’}?  - 
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P (-cos  y)-P  (cos  y))  cos  y 


, ..  " 2 ] r 

+ [P  (-cos  y)  + P^(cos  y >)  sin  y | _ Z B^  cos 


N 


+ P (-cos  y)  - P (cos  Y))sin  Y l ou  sin  w , 


i = 1 


2 . 


/,  , . 4a  sin  vti  (£-fl ri) 

\ *■  ) = 


(P,  (-cos  y)+  P (cos  y))n  cos  y 


" ..  t ^ 2 ^lJ0 

+ ip  (-cos  y)  - P j ( cos  y))  sin  y 1 sin  w.^+ P'  (1)  • n 


2 W, 


i=l 


+ S fp  (-cos  Y)  -P.  (cos  Y))cos  Y+  (p.  (-cos  y)+P  (cos  y) 


. 2 

sin  y 


N , N 

!?i  sin  wi;L+  fP  ^ ( - cos  y ) -P,  (cos  y))  sin  Y Z ou  cos  w^, 


where 


(2.131  l^i  = wfPv(oki)  - Pv(pki)  cos  y sin  wk. 


» ★ 


»»  it  ^ 2 

1J^Pu(pki)  + Pv(pki^  sin  Y cos  Y sin  wki 


M H * 2 

M (p,.  (oki ) _P  i (Pki  0 sin  Y cos  Y sin  wk.  cos  w . , 


Bki  = |,(PY{pki,~t>v<0ki)  sin  3 cos  “ki 


ki 


+ ^(pv(pki)  " Pv(pki})  Sin3y  sin2wki  ' 
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AjJk  - M 0 COS  y |pv(-cos  y)  + (cos  y)J 
+ M0sin2y  |pv(-cos  Y)  - Pv(cos  y)|-  4a2sin  vtt  cos 

N - 2 ? 

+ p Pv^pki^  ^cos  Y + sin  y cos  w .) 

i=l  *-L  ki 

i^k 

N 

+ u . I , pv (pki}  sin2Y  cos2y  (1-cos  w.  . )2 

i=l  v ki 

i A 

r?  * * 2 2 

+ ^ ^ Pv^pki^  (cos  Y ~ sin  y cos  w,  .) 

i=l  *-*-  ki 

i/k 

p £ pv(pki}  siri2y  cos2y  (1+cos  w . )2 

i=l  v ki 

i^k 

+ 2u(cos2Y-sin\)  P^<cos2  -sin2y)-  8p  sin2Y  cos2Y 

"2  2 
Py ( cos  y - sin  y)  , 

B^i  li  cos  y Sin  Y Sin  wki{p>ki>  - P^lp*.)] 
" 3 

+ P Pv^Pk'^  s^n  Y cos  Y (1-cos  w.  .)  sin  w . 

ki  ki 

+ M Pv(pki)  sin3y  cos  Y (1+cos  w ) sin  w,  . , 

Ki  ki 

Aki  = ~M  Pv(pki)  (sin2Y-(-  cos2y  cos  w ) 

KI 

” 2 2 

+ M Pv^Pki^  sin  Y cos  Y (1-cos  w )2  - 

KI 
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i * ,2  2 

~ M pv  C P ) (sin  y-  cos  y cos  w^) 


2 2 ? 

~ M Py(Pki)  sin  Y cos  y (1+cos  w .. ) , 


and 


2 

. 2 

II 

•H 

- cos  y 

- sin  y 

cos 

Wki  ' 

* 

2 

. 2 

ki  = 

cos  y 

- sin  y 

cos 

Wki  ' 

C = Bq  cos  aQ  , 


H = 60  sin  aQ  , 


As  we  can  see,  the  linearized  equations  of  motion  have 

constant  coefficients  and  the  usual  exponential  stability 

analysis  is  possible  for  all  cases.  That  is,  the 

quantities  , 6,  , £ , n can  be  sought  as  a sum 

Xft 

of  exponential  functions  of  the  form  e , l = 1, . . . ,2  (N+l)  . 
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2.3.  Linear  Stability  Analysis  of  the  Four  Cases. 

Case  a.  N Vortices  on  a Circle  of  Colatitude, 

No  Boundary  Condition  at  the  Equator. 

If  we  set  Mq  = 0 in  equations  (T..3),  (2.4),  (2.7) 
and  eliminate  equations  (2.5),  (2.6),  which  involve  the 

coordinates  of  the  polar  vortex,  we  have  the  linearized 
equations  for  the  present  case. 

These  equations  can  be  written  in  the  form 


a ' 

A B >1 

• a ' 

■ a ' 

= 

= L 

k 6 - 

>■  C A ‘ 

3 ^ 

^ 3 - 

where 


and 


a = 


a. 


N 


6 = 


A " (aki}  ‘ [ . 2 


- Ra 
Bki 


4a  sin  vtt  sin  y 


N 


- 

2 . 

l4a  sin  vtt 

; akk  = ° ' 


(2.14)  B = (bki)  = 


A?. 

ki 


2 

4a  sin  vtt  sin  y 


C = (cki.)  = 


B 


ki 


4a^sin  vtt 


N 

' ckk  = " Cki 


i=l 

i/k 


where  the  matrices  A,  B,  C are  of  order  N. 

Furthermore  A,  B,  C are  circular  matrices,  so  that 
it  can  be  shown  that  [8]  the  eigenvalues  of  L are  given  by 


(2.15a) 


X . = a . + (b  e . ) 
3 3-33 


1/2 
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where 


(2.15b) 


a . 

3 


b . 
3 


r (k-l)j 

an  + J2  aik  “ ' 3 


— 1/.../N; 


bU  + j2  bjk  “ 


N 

C-i  = C11  + 2 C"’-  W 

3 1 k=2 


'jk 


(k-1) j 
(k— 1 ) j 


u = 


2i?r/N 
e , 


and  aj'bj'cj  are  the  eigenvalues  of  A,B,C. 

From  (2.14)  we  have  the  following  properties  of  aik 'bik * cik‘ 


an  = 0 

for  all 

values  of  N 

If 

N is  even 

al,N/2  +1  ' 0 ' 

a12 

alN  ' * • * ,al,N/2  ‘ al 

,N/2  +2 

b12 

blN  ' * * * ,bl,N/2  = bl 

,N/2  +2 

C12  = 

C1N  • • * ‘ • c1,n/2  = C1 

,N/2  +2  * 

If 

N is  odd 

0) 

to 

II 

1 

alN  ' * ‘ * ' al,  [N/2 ] +1  = 

"ai,[N/2]+2  ' 

II 

CM 

f—i 

-Q 

blN  '***'  bl,  [N/2 ] +1  = 

bl,  [N/2] +2  ' 

C12 

C1N  '**•'  °1, [N/2 ] +1  = 

Cl,  [N/2 ] +2  * 

Therefore,  we 

get  for 

a^,bj,Cj  the  following  expressions. 
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a . 
3 


= 2i  l 


k=2 


alk  Sin  (wkl  ‘ j) 


(2.16)  b. 

3 


bn  + 2 J b^k  cos  (^^‘j)  + a(N)  cos  7T j b 


k=2 


c . 
3 


4. 

2 I clk(cos  - 1)  + a(N)(cos  7T j —1 ) c , 

k= 2 1 < 


where 


“ki  - (k-1)  it 


f 


P = 


N 

2 

N+l 


for  N even 
for  N odd 


a(N)  = 


1  for  N even 
, 0 otherwise 

Using  (2.13),  (2.14)  and  (2.16)  we  get 


2  . p 

(2.17)  4a  sin  vtt  a.  = -2ip  cos  y l p'(p,.)  sin  to.. 

3 k=2  ^ *1  kl 


sin  (ukl*j)  - 2ip  cos  y sin2y  £ p"  (p  ) sin  to.  . (1- 


k=2 


v 'Kkl‘ 


kl 


•sin  (io 


(2.18)  4a2  sin  vti  b . = - 4a2  sin  vtt  — os  .K  q 

3 sin  y 


kl 


+ 2M  k=2  P^(Pkl)  \^~y  (1_cos  ajkrcos(u)ki*  j)) 


+ sin  yfcos  to.  , -cos  (to.  *j) 


,)} 


kl  VJJkl 

2 p 

+ 2P  sin  y cos  y [ p"  (p  ) ( 1-cos  to.,)2(l+cos  (to,  -j 

k=2  v k-L  kl 


+ a(N)  P‘jPv(Pn/2  +lfl)  " sin  y)  + 


, [N/2  ] +1 
[N/2J+1 


COS  U)^) 

• j ) . 
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" 2 ) 

+ 4 Pv<Pn/2  +1  Sin  Y COS  Y?  (1  + COS  7T  j ) 


2 P 
(2.19)  4a  sin  vttc.  = 2psin  y £ P (p  , ) (cos  (w,  . • j ) -1)  cos 

3 k_2  v k J.  Ki 


w 


kl 


3 r " 2 

+ 2p  sin  y l pv(Pkl)  sin  wkl (cos (wki* j > - 1) 
k 2 


- a (N)  p P (p 


v n/2 


+1  sin  y (cos  nj  - 1) 


(2.20)  fi  = 


2p 

. 2 . 

4a  sin  vtt 


r P 


[ Py(Pk  cos  Y (1-cos  wkl) 
k — 2 


♦ o(N)  cos  Y Pv(0n/2  ♦ !.  1> 


A study  of  the  sign  of  c.,  shows  that  c.,  < 0 for 

IX  IK 

all  values  of  y,  N,  k,  cp  whenever  p > 0 , which  implies 

0^  >_  0 for  all  j,  y » N,  cp. 

Therefore,  we  have  from  (2.15a,b)  that  if  b^  ^ 0 

for  all  j,  L has  only  purely  imaginary  eigenvalues,  and 

if  bj  > 0 for  some  j,  then  L has  complex  eigenvalues 

(except  when  c^  = 0) , the  second  alternative  would  mean 

exponential  instability  for  the  linearized  equations  if 

the  real  part  of  at  least  one  eigenvalue  is  positive. 

p can  be  factored  from  a.,  b.,  c.;  let 

3 3 3 

aj  = "aj  ' bj  = ' c:  = ’S 

then 


* “ 


* .it 

1 + /bZcZ 


'Z  - 
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therefore  p affects  only  the  magnitude  of  A and  does 
not  affect  in  any  way  the  stability  properties  of  a given 
configuration  of  vortices. 

If  P < 0 we  have  to  reverse  the  inequalities  of  the 
last  paragraph. 

We  can  state  now  the  linear  stability  analysis  reduces 

to  finding  those  values  of  N,  0,  c that  will  yield  a 

P 

neqative  value  of  b_^  when  p > 0 or  a positive  value 
of  b^  when  p < 0 for  all  values  of  j . Such  values  of 
W,  0,  c will  define  the  linearly  stable  range.  From 
(2.15a) , we  see  that  each  value  of  j gives  two  eigenvalues 


of  L. 

If  j = N,  we  have  a^  = c^  = 0,  which  means  that  f>  is 
an  eigenvalue  of  L with  multiplicity  2.  Summing  eq.  (2.3) 
and  (2.4)  with  pg  = 0 , we  have 


Z 3 

k=  1 


N N 


N N 


k * l Z aikBi  + Z l cik(ctk"ai) 

K k=l  i=l  lk  1 k=l  i=l  lk  k 1 


iA 


iA 


N r N 


[ Z a,,.  |B,  = 0 


Integrating  we  get 


i=l  M<=1 
Mi 


N 


ik  i 


Also 


k=  1 
N N 


Z B.  = constant 


N . N N N N 

Z \ = Z Z aik(A'V  + Z Z bikek 

:=1  K k=l  i=  1 1K  k 1 k=l  i=l  lK  k 


N N 


N f N 


= l I bi  ? - l l b 

k= 1 i= 1 11  k k=l  i= 1 


V 


N 
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H r — i S 


3y  an  appropriate  choice  of  displacements  £ 

K 

we  can  obtain 


(2.21) 
and  then 
(2.22) 


N 


l B 

k=l 


k 


0 


c = constant  . 


These  invariants  can  also  be  obtained  by  linearizing 
and  f2  respectively  defined  in  (1.14),  (1.15). 

If  we  make  j = 1 in  (2.17),  (2.18),  (2.19)  we  have 


(2.23)  4a2sin  vtt  a 


. 2 


= -2i v cos  Y [ P (p.,)  sin 


U) 


k=2 


kl 


. 2 P 


. 2 


- 2i U COS  Y sin  Y [ P (PRl)  sin  wvl  (1-cos  u>ul  ) 


k=2 


kl 


kl' 


(2.24)  4a  sin  vtr  b^  = 


-4a^sin  vtt 

sm  y 


+ 2p  j2  p>kl>  IlHT  si"2  “kl 


+ 2p  sin  y cos  Y / P.,(Pkl)  sin  “vi (1‘cos  ) 


k = 2 


kl 


(2.25)  4a  sin  vtt  c1  = 2u  sin  y [ Pv{pkl)  (cos  aJkl"1)  cos  w 


k-2 


. 3 


k = 2 


v'Mkl‘ 


Jkl 


ukl 


-cos 

“ki> 

kl5 

(cos  aj. 

-1) 

kl 


+ a(N)  2v  p;(pn/2  +1,  1>  sin 
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Using  (2.23)  and  (2.24)  in  (2.25)  we  get 


, 2 . i sin  y ir,  • \ 

4a  sin  vt  c,  = 4a  sm  vtt  — — > (ft-ia..) 

1 cos  y 1 


Using  (2.23)  in  (2.24), 


4a2  sin  vtt  = -4a2sin  vtt  • 


Hence 


X^  = a^  + (ft-  ia^ ) i = 


- ift  and 


2a^-ift 


For  j = N-l,  one  of  the  eigenvalues  is  = i^#  i 

the  period  ft  is  an  eigenvalue  of  L. 

For  y = tt/2  , we  have  for  a fixed  c^  and  N, 


a_.  = 0 for  all  j = 1,...,N 

P 

bj ‘ 2u  J , Vpki>  (cos  “kr003 

- \i  a (N)  P^(Pn/2  +i,  i}  (1  + cos  "3* 

P 

c.  = 2u  l Py(Pkl>  (cos  “»kl*j  _1)  cos  wki 
-1  k=2 

- a(N)  p Pv(pn/2  +1,  1}  ({"1)3  ' L)  * 


It  follows  that 


b,  = 0 

L 

b„  . = 0 which  means  ft  = 0 . 
N —1 
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From  the  above  considerations,  we  can  conclude  that 
when  y = tt/2  , we  have  only  zero  eigenvalues  for 

configurations  with  2 and  3 vortices. 

As  we  go  closer  to  the  pole  ft  increases  and 


00  if  u > 0 

lim  ft  = ■ 

0-0  -«>  if  v < 0 

when  y - 0,  then  p..  - -1  for  all  i,k  = 1,...,N. 
As  the  argument  tends  to  -1, 


Pv(x) 


sin  vtt 

TT 


2 

x+1 


and 


P^(x)  'v 


sin  vtt  2 

(x+1)2 


Using  these  approximations  of  P'(x)  and  P (x)  in 
the  expressions  for  ay  b , c ^ (2.17-19)  the  factor 
sin  v 71  cancels  and  we  can  conclude  that  for  small  colatitudes 
the  eigenvalues  do  not  depend  on  the  coriolis  parameter. 

The  numerical  results  listed  in  the  Appendix  suggest 
the  following  conjecture. 

Given  a configuration  of  N vortices,  there  is  a colatitude 
Y,  that  depends  on  N,  such  that  for  0 <_  y the  configuration 
is  linearly  stable  and  for  0 > Y is  exponentially  unstable. 
This  particular  value  of  y decreases  for  increasing  N. 


Case  b.  N Vortices  on  a Circle  of  Colatitude , 1 Polar  Vortex , 
No  Boundary  Condition  at  the  Equator. 

Equations  (2.3),  (2.4),  (2.5),  (2.6)  give  us  the 
linearized  equations  for  the  case  of  N vortices  on  a circle 
of  colatitude  G and  a polar  vortex  with  no  boundary 
condition  at  G = tt/2  . 

These  equations  can  be  written  as 


The  addition  of  a polar  vortex  destroys  the 
circulant  property  discussed  before  and  the  eigenvalue 
calculation  has  to  be  carried  out  by  an  iterative 
procedure . 

We  have  to  consider  a new  parameter,  namely,  the  polar 
strength  The  quotient  plays  an  important  role 

in  the  linear  stability  analysis. 

Some  of  the  general  properties  which  arc  valid  for  all 
values  of  the  parameters  are  the  following: 
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(a)  the  angular  velocity  £i,  as  given  by  (1.7)  is  an 
eigenvalue  of  L • 

It 

(b)  L has  a double  zero  eigenvalue. 

(c)  lim  H ‘J0  V0)  * ° ' 

Y-*tt/2  4a  sin  vn 

In  the  derivation  of  the  linearized  equations  (2.3) -(2.6) 
mixed  spherical  and  polar  coordinates  wer-  used  to  describe 
the  motion  of  the  polar  vortex  and  the  vortices  on 
the  circle  of  colatitude.  This  was  necessary  in  order 
to  give  coordinates  to  the  polar  vortex  so  that  the  equation 

of  motion  is  lonsingular  at  M = °* 

If  we  sum  equations  (2.3)  and  (2.4)  over  k or 

linearize  ( 2 . 14  ) , ( 2 . 15 ) , we  get 


? 6,  = constant  and  I *k  = constant. 

k-1  " k=1 

That  is,  we  have  two  linear  invariants  of  the  motion 
but  they  are  independent  of  the  position  of  the  polar  vortex 
From  equation  (1.7)  we  have  that  if 


N 

cos  Y l p'  (Pij)  (l"coS  "li 
UQ  i=2  * 

0 P (-cos  Y) 

V 

then  fl  = 0.  Since  fi  is  an  eigenvalue 
we  have  a quadruple  0 eigenvalue  for  all 

If  U > o , u0/u  1 C%  n 1 0 and  the 

positive  direction  (West  to  East)  . 

n < 0. 


) 

* 

- = c 

of  Lp,  for  U0/U  = 
values  of  N,  Cp,  0 
rotation  is  in  the 
If  u > 0,  u0/w  < 


* 

C 
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Some  results  on  the  algebraic  eigenvalue  problem 


for  N - 3 is  presented  in  this  section.  For  other  values 

of  N,  corresponding  results  are  given  in  the  Appendix. 

The  calculations  were  done  using  u = .15  • 10  m /sec.; 

values  of  c = .2,  ...,  1.0;  0 = 5° ,6° , . . . , 90°  and 

Uq/u  < 0.  They  can  be  summarized  as  follows. 

N = 3.  uQ/u  < 0.  For  a(cp)  < PQ/p  < 0, 

max  'a(c  ) | - 2.0  we  have  stability  in  the  entire  Northern 
c p 

hemisphere.  For  | m q /m | > a(c  ),  the  stability  region  is 

reduced  to  a cap  that  decreases  with  increasing  | y /\i  \ . 

In  the  table  below  we  have  the  limit  0 of  the  cap  for 

c 

different  values  of  Pq/u* 


Table.  Of  0 for  N = 3 
c 

Configuration  is  stable  for  0 <_  0 0c> 


For  Pq/u  = C we  have  linear  stability  for  all  0,  c 
U /u  > 0.  For  0 < u /m  < b ( c ) , max  b(c  ) *•  0.75,  we 

U v P c r* 

have  linear  stability  for  0 < 9 <_  n/2.P 

If  Vq/v  increases,  a belt  of  instability  appears,  a: 
the  width  of  this  belt  increases  with  increasing  uQ/\i. 

In  the  table  given  below,  we  have  the  limits  of  the 
instability  belt,  0.  and  °2  , unstable  for  01  < 6 ' 02 
stable  otherwise,  for  some  values  of  c and  Wq/)j* 


* 

Table 

Configuration  is  unstable  for  l)  < 9 < 92. 


1 

~0~1 

1 

25 

1 

5 

2 

0 

2 

5 

3. 

0 

4 

0 

5 

0 

C N. 

p X, 

b 

°2 

hi 

°2 

r* 

1 

n 

2 

A 

1 

•H2 

\\ 

2 

h 

_^2_ 

h 

il 

. 2 

0 

0 

11 

28 

9 

37 

9 

47 

8 

52 

8 

56 

7 

60 

7 

63 

.3 

0 

0 

10 

29 

9 

39 

9 

49 

8 

54 

8 

58 

7 

62 

7 

65 

.4 

0 

0 

10 

32 

9 

41 

8 

51 

8 

57 

8 

60 

7 

6 4 

7 

67 

.5 

0 

0 

10 

34 

9 

44 

8 

5 4 

8 

60 

8 

63 

7 

67 

7 

69 

.6 

0 

0 

10 

37 

9 

46 

8 

57 

8 

62 

8 

66 

7 

70 

7 

72 

. 7 

14 

17 

10 

39 

9 

49 

8 

59 

8 

65 

8 

68 

7 

72 

7 

74 

. 8 

13 

22 

10 

42 

9 

51 

8 

62 

8 

67 

8 

71 

7 

75 

7 

77 

.9 

12 

25 

10 

43 

9 

53 

8 

64 

8 

69 

7 

73 

7 

77 

7 

79 

1 .o 

12 

28 

10 

44 

11 

55 

8 

65 

8 

71 

7 

; 

75 

7 

79 

7 

81 

♦Entries  in  table  are  in  degrees. 


Using  these  tables,  we  can  construct  the  curves 
of  linear  stability  in  the  plane  (0,6q/6)  for  a fixed  c^ 
See  Figure  A-70  in  the  Appendix. 


Case  c.  N Vortices  on  a Circle  of  Colatitude, 


Boundary  Condition  at  the  Equator. 


Setting  uQ  = 0 in  equations  (2.9),  (2.10),  (2.13)  and 
eliminating  equations  (2.11),  (2.12)  we  have  the  linearized 
equations  for  the  configuration  to  be  considered  now. 

The  equations  can  be  written  in  the  form 


f a 1 

f A2 

B2 

r a 

l 0 > 

l C 

Ao 

l 6 

= lb 

[ J 

and  using  (2.9),  (2.10)  with  Uq  = 0 we  have  that 


A2  * (aki>  - 


-B 


a 

ki 


l4a  sin  vir  sin  y 


( ~ A 


ki 


l4a  sin  vtt 


' akk=0 


(2.26)  B2  = (bRi)  = 


C2  - «cki>  * 


a“ 

ki 


i.2. 

v4a  sin  vtt  sin  y 


B» 

ki 


. 2 . 

4a  sin  vtt 


N 


c,.,.  - - l 

i=l 

iA 


' "kk  " A,  Cki  ' 


A2'B2'C2  are  c;*-rculant  matrices  and  therefore 
equations  (2.15)  also  hold  here. 

Using  (2.16),  we  can  get  equations  similar  to  (2.17)— 

(2. 20)  . As  before,  we  can  observe  that  c^  >_  0 when  y ^ 0 
and  conclude  that  if  b^  > 0 for  some  j,  we  have  exponential 
instability  and  if  b^  <_  0 for  all  j,  we  have  a linearly 
stable  configuration.  In  this  case  0 is  also  an  eigenvalue 
with  multiplicity  two. 
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Invariant  expressions  similar  to  (2.21)  and  (2.22)  can 
be  obtained  for  the  configuration  being  discussed  now. 

The  behavior  of  a ^ , b_.,  c^  as  y *♦  0 is  similar  to  that 
of  case  a and  we  can  conclude  also  that  given  N there  is 
a y (N)  such  that  for  y <_  y (U)  we  have  linear  stability. 

For  small  colatitude,  the  eigenvalues  do  not  depend 
on  the  coriolis  parameter;  this  can  be  seen  by  using  the 
argument  given  in  case  a. 

Even  though  the  above  properties  are  similar  to 
cases  a and  b,  the  addition  of  a boundary  condition  at 
the  equator  causes  some  fundamental  differences. 

First  of  all,  ft  as  given  by  (1.12)  with  p^  = 0 
is  not  an  eigenvalue  of  . We  have  that 


ft  = 


; 

4a  sin  vtt 


p cos  y 


N 

I 

i = 1 

i/k 


Pv(pki)  (1-COS  wki) 


(1+cos  w..)  + 2p'(cos  2y) 

V 


as  y tt/2  , the  terms  under  the  summation  sign  tend  to  0 
because  of  the  factor  cos  y.  But, 

lim  cos  y p'(cos  2y)  = 

Y+tt/2 

so  that 

lim  ft  = * if  p > 0 
Y-*-tt/2 

and 

lim  ft  = -oo  if  p < 0 . 

Y +tt/2 
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On  the  other  hand 

lim  fi  = ± 00  for  p < 0 
Y+0 

Secondly,  as  we  will  see  later,  the  zones  of  stability  in 
the  plane  (c^,0)  are  reduced  by  the  existence  of  a 
boundary  condition. 

Results  from  the  numerical  calculations  are  summarized 
in  the  Appendix. 
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Case  d.  N Vortices  on  a Circle  of  Colatitude , One  Polar 


Vortex,  Boundary  Condition  at  the  Equator. 

Equations  (2.9)  - (2.12)  are  the  linearized  equations 
of  motion  for  the  case  to  be  discussed,  i.e.  (N+l)  vortices 
satisfying  the  boundary  condition  v = 0 at  0 = tt/2. 

The  equations  can  be  expressed: 


• 

5 ^ 

E3 

°3 

6 

6 • 

• 

a 

■ 6 ■ 

F3 

k 

A3 
C 3 

B3 

A3 

a 

6 ■ 

1 

= L 

P 

a 

■ B • 

The  matrix  L is  not  circular  and  the  calculation  ot 
P 

the  eigenvalues  has  to  be  carried  out  by  one  iterative 
procedure . 

As  in  case  b,  we  have  a new  parameter  Uq  and  the 
linear  stability  properties  depend  on  the  value  of  yQ/y. 
Some  of  the  general  properties  are  the  following. 

I 

(a)  the  angular  velocity  ft  is  not  an  eigenvalue  of  L^. 

(b)  Lp  has  a double  zero  eigenvalue. 

(c)  lim  ft  = 00  if  y > 0 and  lim  ft  = -00  if  y < 0. 

Y+it/2  Y+tt/2 

The  invariants  are  also  independent  of  the  position 

of  the  polar  vortex  and  mixed  coordinates  were  used  to 

describe  the  linearized  motion  of  the  polar  vortex  and 

the  vortices  on  the  circle  of  colatitude.  From  equation 

(1.12),  we  see  that  if 
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★ 


u 


P(-cos  Y)  + P 1 ( cos  Y) 


= c 


/ 


then  . = 0 and  the  vortices  remain  stationary. 

* # 

For  Mq/u  ± c and  UQ/u  < c the  rotation  is  West 
to  East  and  East  to  West  respectively.  Since  ft 
is  not  an  eigenvalue  of  i/  , we  do  not  have  a quadruple 
zero  eigenvalue  for  uQ/u  = c*  as  in  case  b. 

As  in  the  previous  case  with  polar  vortex,  some  results 

are  given  for  N = 3.  In  the  Appendix,  results  for  other 
values  of  N are  given. 


N = 3-  VQ/U  < 0.  For  negative  values  of  uQ/u  we  have 
a cap  of  stability.  Given  below  are  values  of  the  limit, 
°c  ' this  cap.  (See  the  table.) 
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Table 


Values  of  6 . 

c 

Configuration  is  stable  for  0 <_  0 £ 0 . 


0.2 

-.5 

-1.0 

t 

-1.5 

■ 

-2.5 

m 

-3.0 

D 

42° 

39° 

35° 

30° 

23° 

39° 

n 

43° 

42° 

39° 

35° 

30° 

23° 

39° 

41° 

43° 

42° 

39° 

35° 

31° 

24° 

40° 

37° 

43° 

42° 

39° 

36° 

31° 

24° 

41° 

35° 

44° 

42° 

40° 

36° 

32° 

24° 

42° 

34° 

43° 

40° 

, 37° 

32° 

25° 

n 

34° 

40° 

37° 

32° 

25° 

B 

34° 

41° 

37° 

33° 

26° 

35° 

y 

I 

41° 

1 . 

38° 

1 

34° 

26° 

44° 

| 

36° 

Table . 


The  belt  disappears  for  increasing  values  of  Uq/v  and  for 
vJg/u  = 2.0,  2.5,  we  have  again  a cap  of  stability. 


Configuration  is  stable  for  0 < 0 < 0 . Values  of  0 : 

- - c c 


V 

Increasing  Mq/p  > the  new  cap,  is  broken  into 
another  cap  and  a belt  of  stability.  Given  below  are 
some  values  of  the  limit  0 of  the  cap  and  0^,02  the 
belt  for  different  values  of  Pq/m* 

Configuration  is  stable  for  O<0<G  , 6,  <0  <0... 

— — c 1 — — 2 


V“ 

3.0 

4.0 

6.0 

cap 

belt 

cap 

c 

E 

0 

c 

V 

02 

0 

c 

°1 

°2 

.2 

7° 

SB 

63° 

6° 

63° 

67  0 

6° 

m 

73° 

. 3 

7° 

63° 

63° 

6° 

63° 

67  0 

6° 

67° 

73° 

.4 

7° 

64° 

64° 

6° 

62° 

68° 

6° 

O 

VO 

VO 

73° 

.5 

7° 

- 

- 

6° 

65° 

68° 

6° 

66° 

73° 

.6 

7° 

- 

- 

6° 

66° 

68° 

6° 

65° 

72° 

.7 

7 

- 

- 

6° 

67° 

68° 

6° 

65° 

72° 

. 8 

7° 

- 

- 

6° 

67° 

6° 

67° 

72° 

.9 

7° 

- 

- 

6° 

67° 

O 

00 

VO 

6° 

68° 

72° 

1.0 

7° 

- 

- 

6° 

68° 

Oh 

00 

o 

6° 

Oh 

00 

o 

71° 

If  we  compare  these  tables  with  those  of  case  b,  it 
is  clear  that  the  boundary  condition  reduces  the  zones  of 
exponential  stability  in  the  plane  (6,Ug/u). 

The  curves  of  linear  stability  for  N = 3 and  other 
values  of  N are  given  in  Figures  A-69a,b,c  of  the  Appendix. 
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3 . Numerical  Solution  of  the  Nonlinear  Equations  of  Motion 


3. 1 Introduction 

The  linear  analysis  of  Section  2,  established  certain 
zones  of  exponential  stability  and  instability  in  the 
( 9 , Uq/u)  plane  and  (0,cp)  plane  for  configurations  with 
and  without  polar  vortex  respectively.  But,  for  the  entire 
range  of  parameters,  zero  is  an  eigenvalue  of  multiplicity 
two  and  in  some  cases  of  multiplicity  four. 

It  is  therefore  desirable  to  study  the  type  of 
solution  associated  with  the  multiple  zero  root,  and  its 
effect  on  the  linear  stability  properties,  before  describing 
the  results  of  the  numerical  integration  of  the  nonlinear 
equations  of  vortex  motion. 

First  let  us  analyze  configurations  without  a polar 
vortex.  As  in  a previous  section,  we  linearize  about  an 
equilibrium  solution: 


(3.1) 


0k<t)  = Y = constant  , 

*k(t>  = (k— 1 ) ^ + at  , 

K N k = 1,2, 


,N. 


That  is,  a perturbation  of  the  solution  (3.1)  is  considered 
to  have  the  form 


(3.2) 


9k(t)  = Y + 6k(t)  » 

♦k(t)  = (k-1)  ^ + nt  + cxk(t)  , 

k — 1,*..,  N. 


where  ak  and  0^  are  small  at  time  t = 0 . 
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If  the  solution  (3.2)  remains  close  to  the  equilibrium 
solution  (3.1)  for  all  t,  that  is,  if  ot^tt)  and  8^(t), 
k = 1,...,N,  remain  small  for  all  t,  we  say  that  the 
conf iquration  is  stable. 

Assume  that  all  the  eigenvalues  of  the  linearized 
system  (2.14),  (2.26)  have  negative  or  zero  real  part,  i.e., 
the  configuration  is  exponentially  stable.  By  using  the 
invariants  defined  in  Section  2, 


N 

(3.3) 

I 

k=l 

6k(t) 

= c^  = constant 

N 

• 

(3.4) 

l 

k=l 

“k(t) 

= c^*c2  = constant  , 

we  see 

that  if  c^ 

= 0, 

then  the  solutions  a^(t),  Bk(t), 

k = 1,...,N,  are  bounded  for  all  t,  that  is,  the  time 
increasing  components  of  the  form  (constant) *t  associated 
with  the  double  zero  eigenvalue  are  eliminated  by  the 
condition  c^  = 0 . 

If  c^  ^ 0,  an  integration  of  (3-4)  gives 
N 

“k(t)  “ C1  c2  ‘ + c3  ' 

so  that  if  c^  5*  0 , then  the  multiple  zero  root  yields 
time  increasing  component  in  the  longitudinal  direction. 
According  to  the  previous  remarks  about  stability,  such 
a linear  growth  would  have  to  be  termed  an  algebraic 
instability . 
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N 

In  Section  2,  we  saw  that  c2  = [ b.^  , where  b^  is 

i = l 

qiven  by  (2.14).  By  using  (2.7),  (2.14),  and  (1.7),  it  is 
easy  to  prove  that 

- y b = — 
c2  " . - . bii  * 

i=l 

Hence,  the  change  of  the  longitudinal  components  is 
proportional  to  the  variation  of  the  angular  velocity  with 
colatitude.  This  is  a reflection  of  the  fact  that  when  the 
vortices  are  slightly  displaced  in  colatitude  initially, 
the  angular  velocity  of  the  equilibrium  solution  changes. 

Note  that  the  change  in  angular  velocity  is  of  the  order  of 

3Q 

cr  3y  * 

Furthermore,  as  we  will  see  later  in  this  section  and  in 
the  Appendix,  whenever  c^  f 0,  the  vortices  tend  to 
oscillate  about  a new  equilibrium  position,  having  a slightly 
modified  angular  velocity  of  rotation. 

It  seems  reasonable  after  the  considerations  just 
discussed  to  define  stability  by  saying  that  the  configura- 
tion is  stable  if 

S (t) , k = 1 , . . . ,N , remains  small  for  all  t; 
k 

and  if  ak(t),  k = 1,...,N,  remains  small  for  all  t; 
after  subtracting  a term  of  the  form  (constant) • t. 

When  the  configuration  has  a polar  vortex,  the 
expressions  (3.3),  (3.4)  are  independent  of  the  coordinates 

of  the  polar  vortex.  In  this  case,  the  effect  of  the  zero 
roots  on  the  solution  is  not  clearly  defined. 
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Nevertheless,  if  only  the  circle  vortices  are  perturbed 
initially,  then  the  nonlinear  calculations  show  that  the 
preceding  qualitative  description  of  the  motior  of  the  non- 
polar vortices  also  applies  to  their  behavior  in  this  case. 

On  the  other  hand,  if  only  the  polar  vortex  is  displaced 
initially,  then  the  qualitative  behavior  of  the  resulting 
motion  cannot  be  so  easily  described. 

The  numerical  solutions  of  the  nonlinear  initial  value 
problem,  will  give  us  an  overall  picture  of  the  motions  in 
the  regions  of  linear  stability  and  instability  and  most 
important  of  all,  it  will  give  some  indication  of  the 
validity  of  the  linear  analysis. 

The  numerical  integrations , were  carried  out  for  a 
period  of  4000  hours  of  physical  time. 

Several  types  of  initial  perturbations  were  considered 
in  the  numerical  calculations: 

(a)  Displacement  of  the  polar  vottex  in  the  6 direction. 

(b)  Displacements  of  one  circle  vortex  in  the  positive 
and  in  the  negative  0 direction 

(c)  Disr lacements  of  the  circle  vortex  in  the  positive 
and  in  the  negative  <J>  direction. 

The  disturbance  produced  by  displacing  only  the  polar 
vortex  is  more  nearly  symmetrical  than  is  a disturbance 
produced  by  an  initial  displacement  ir  the  <J>  direction  of 
a circle  vortex.  In  fact  the  magnitude  of  the  latter 
disturbance  is  of  the  order  of  50%  larger  than  the  magnitude 
of  the  former  disturbance,  for  equal  initial  angular  displacement. 
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In  general,  it  is  possible  to  classify  the  motions 
obtained  by  numerical  integration,  as  follows: 

1.  Stable  linear  motions.  This  type  of  motion  is 
typical  of  con figuration  with  parameter  values  in  the 
linear  stable  zone.  They  are  characterized  by  the 
fact  that  each  component  (0k,$k)  of  the  motion  is 
bounded,  at  all  times,  by  the  initial  perturbation 
i.e.,  the  amplitudes  of  the  oscillations  are 
proportional  to  the  amplitude  of  the  initial  perturba- 
tion (a  linear  scaling) . 

2.  Nonlinear  periodic  motions,  followed  by  an  unstable 
departure  away  from  the  initial  equilibrium  configura- 
tion. These  types  of  motions  occur  mostly  for  parameter 
values  in  the  exponentially  unstable  zone.  The  nonlinear 
periodic  motion  persists  for  a finite  interval  of  time 
that  depends  mainly  on  u /u  and  c . 

The  amplitudes  of  the  subsequent  oscillations  are  very 
large,  as  compared  with  the  initial  perturbation,  and  they 
are  not  related  in  a linear  way  to  the  size  of  the  initial 
disturbance. 

Configurations  consisting  of  3 circle  vortices  and  one 
polar  vortex,  with  and  without  boundary  conditions  are 
treated  in  this  section.  Configurations  without  a polar 
vortex  are  studied  in  the  Appendix. 
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3 . 2 3 Vortices  of  Equal  Strength  y , Equally  Spaced  on 

a Circle  of  Colatitude  0 , One  Polar  Vortex  with 


Strength  y^  , No  Boundary  Condition  at  the  Equator. 

Here  we  analyze  the  numerical  behavior  of  one  case 
in  detail.  The  stationary  equilibrium  configuration  is 
defined  by  the  parameter  values 

0 = 60°,  c = 0.2,  yA/y  = -.597963...  . 

P 0 

According  to  the  results  given  by  the  linear  analysis, 
this  configuration  is  in  the  exponentially  stable  region, 
and  zero  is  an  eigenvalue  of  multiplicity  four. 

In  Figures  3.3  and  3.2,  the  components  of  the  motion 
are  shown  for  a polar  perturbation  of  the  form 
(e60  = 5.0°,  eotg  = 0.°)  and  (e3q=  3.0°,  eaQ  = 180°) 
respectively. 

In  Figure  3.1,  we  observe  a steady  increase  of  the 
polar  components  in  time.  This  is  characteristic  of  the 
motions  that  result  from  initial  displacements  of  the  form 
cBq  ^ 3°,  coiq  = 0°.  On  the  other  hand.  Figure  3.2  shows 
that  a stable  motion  with  recurrent  trajectories  occurs  for 
an  initial  displacement  in  the  direction  cBq  = 3°,  eotQ  = 180°. 
This  asymmetry  of  the  response  for  displacements  in  the 
eaQ  = 0°  and  eaQ  = 180°  directions  arises  from  the  asymmetry 
of  the  displacements  relative  to  the  positions  of  the  circle 
vortices  . 

Figures  3.3  and  3.4  show  the  resulting  motions  for  an 
initial  displacement  of  the  circle  vortex  at  <j>  = 0°,  in  the 
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positive  and  negative  0 directions  respectively.  In  both 
cases  we  find  that  the  longitudinal  components  increase 
linearly  in  time.  For  small  initial  displacements,  say, 
-10°  <_  e8^  <_  10°,  the  change  in  the  longitudinal  components 
predicted  by  the  linear  theory  (see  (3.1))  is 


c 


3ft  . 
1 3y 


6.795  • 10 


dgg 

hour 


] 


the  computed  nonlinear  change  is  observed  to  be  of  the  order  of 

e0.  • 7.5075  • 10‘4  [-9^  3 . 

1 hour 


Clearly,  the  nonlinear  effects  are  small.  In  this  case, 
the  linear  theory  gives  a good  qualitative  description  of 
the  motion  for  !^8^|  <_  10°. 

In  Figure  3.5,  we  have  an  example  of  the  response  of 
the  configuration  to  an  initial  displacement  in  the  negative 
<P  direction;  that  is  ea^  = -2°.  Observe  that  there  is  no 
steady  drift  in  the  longitudinal  direction,  but  there  is 
some  evidence  of  increasing  components  at  the  pole.  The 
majority  of  the  integrations  performed  show  that  an  initial 
perturbation  of  a circle  vortex  in  the  longitudinal  direction 
produces  a "greater  deviation"  from  the  equilibrium  state 
than  does  any  other  perturbation  of  equal  initial  angular 
displacement. 
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List  of  Figures 


Each  of  Figures  3.1  - 3.5  consists  of  seven  sub- 
figures labelled  a,b,c,...,g.  The  subfigures  are  plots 
of  the  following  quantities. 

a.  Stereographic  projection,  on  the  tangent  plane  at 
the  pole,  of  the  trajectories  of  the  centev"'  of  the 
vortices.  The  trajectories  are  plotted  in  two 
differently  amplified  scales.  One  scale  for  the  polar 
vortex,  the  other  for  the  circle  vortices.  The 
amplifi cation  factors  may  be  found  by  observing 
Figures  b, . . . ,g. 

b.  eQ(t) 

c.  S^t) 

d.  02(t) 

e.  9 3 ( t ) 

f.  (k-1)  ^ - 4>k  (t)  , k = 1,2,3 

g.  <t>0(t) 

For  example,  Figure  3.2b  represents  the  component  0g(t) 
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Figure  No. 

Colatitude  of  the 
Circle  Vortices 

Coriolis 

Parameter 

Initial 

Displacement 

3.  1 

60° 

.2 

caQ=  5.0°,  eaQ=  0° 

3.2 

60° 

.2 

g60=  3.0°,  ea0=  180° 

3.3 

60° 

.2 

c60=  0.01°, 

eaQ=  0.0°,  eB.  = 2° 

3.4 

60° 

.2 

g Bq=  0.01°, 

caQ=  0°,  c81=  -2.0° 

3.5 

Os 

O 

o 

.2 

cSQ=  0.01°, 

EaQ=  0°,  Eai  = -2.0° 
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3 . 3 3 Vortices  of  Strength  \i , Equally  spaced  on  a Circle 

of  Colatitude  0,  One  Polar  Vortex  with  Strength  VJg  , 
Boundary  Condition  at  the  Equator. 

Here  we  study  two  cases  in  detail,  for  the  coriolis 
parameter  cp  = 0.4. 

The  linear  analysis  establishes  that  if  the  ratio 
VJg/y  is  taken  so  that  the  circle  vortices  remain  stationary 
zero  is  an  eigenvalue  of  multiplicity  two.  Furthermore, 
for  0 < 40°,  the  configuration  is  exponentially  stable; 
whereas  it  is  exponentially  unstable  for  0 >_  41°. 

That  is,  numerical  experiments  were  conducted  in  the 

two  cases: 

Case  (a)  6 = 35°,  cp  = .4,  \iQ/\i  = -0.7813...,  i.e., 

in  the  exponentiall  stable  range,  and 

Case  (b)  6 = 45°,  cp  = .4,  Ug/U  = -0.83107...,  i.e., 
in  the  exponentially  unstable  range. 

Case  a.  In  Figures  3.6  and  3.7,  we  have  typical  examples 
of  the  resulting  motions  for  an  initial  polar  displacement, 
that  is,  they  respectively  plot  the  results  for  the  initial 
displacement  (eBg  = 2.0°,  eaQ  = 0°)  and  (eBg=  1.0°,  eotg-  180  )• 
We  found  that  the  motions  were  linearly  stable  motions 

for  initial  displacements  | e 6 Q | 1 5°. 

In  Figure  3.8  we  plot  the  results  for  an  initial  circle 

vortex  displacement. 

When  the  perturbation  is  in  the  positive  and  in  the 
negative  0 direction  (Figure  3.8),  we  note  that  in  this  case 


no  drift  occurs . 


This  phenomenon  can  be  explained  by  the  fact  that  for 


the  configuration  studied 


3f2  -7  1 

- 0.1637  • 10  ' [~±~] 
°Y  sec 


and 


ci  §7-  • “-5  1-^1- 


Therefore,  after  4000  hours  of  numerical  integration, 
the  linear  increase  is  too  small  to  be  noticeable. 

Figures  3.9  and  3.10  show  typical  motions  when  the 
initial  displacements  are  in  the  <J>  direction.  The  ampli- 
tude of  the  perturbation  is  larger  in  the  case  of  a boundary 
condition  than  in  the  case  of  no  boundary  condition. 


Case  (b) . Figures  3.11  and  3.12  describe  the  behavior  of 
the  components  of  the  motion  when  the  initial  configuration 
is  in  the  exponentially  unstable  range.  Even  though  the 
sjflplitude  of  the  initial  displacement  is  small,  say 


e8g  or  - 0.1°,  we  find  that  the  motions  are  nearly 


stable  motions  for  about  1200  hours.  After  that,  we  have 
a rapid  departure  away  from  the  equilibrium  position. 

As  was  observed  in  the  calculations  of  Section  3.2, 
an  initial  longitudinal  displacement  of  a circle  vortex 
results  in  a stronger  perturbation  than  that  which  arises 
from  an  initial  displacement  in  the  latitudinal  direction. 

It  is  interesting  to  note  that  a configuration 
satisfying  a boundary  condition  is  more  sensitive  to 
initial  displacements  than  is  a configuration  that  does 
not  satisfy  a boundary  condition. 
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List  of  Figures 


Each  of  the  Figures  3.6  - 3.12  consists  of  seven 

subfigures  labelled  a,b ,g.  The  subfigures  are  plots 

of  the  following  quantities: 

a.  Stereographic  projection,  on  the  tangent  plane  at 
the  pole,  of  the  trajectories  of  the  center  of  the 
vortices.  The  trajectories  are  plotted  in  two  differ- 
ently amplified  scales.  One  scale  for  the  polar 
vortex,  the  other  for  the  circle  vortices.  The 
amplification  facto  is  may  be  found  by  observing 
Figures  b , . . . , g . 

b.  Oq  ( t) 

c . 0 ^ ( t ) 

d . C>2  ( t) 

e . 0 3 ( t) 

f . (k  — 1 ) ~ - 4>k  ( t)  , k = 1,2,3 
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4.  Numerical  Method 


The  linearized  eigenvalue  calculation  of  Section  2 
was  carried  out  by  using  the  eigensystcm  subroutine  package 
EISPACK  [9].  The  routines  from  this  package,  reduce  the 
input  matrices  (L,  Lp  , , Lp  of  Section  2)  to  upper 

Hessenberg  form  by  means  of  similarity  transformations  and 
then  computes  the  eigenvalues  by  applying  the  QR  method. 

The  numerical  integration  of  the  nonlinear  equations 
of  vortex  motion  was  carried  out  by  using  the  modified  Euler 
method  [5]  defined  as 


* t t 

u =u  + At  f ( t , u ) 


t+At  t At  r,,  t * . 

u = u + — (f(t,u  ) (-  f ( t+A t , u )} 


The  equations  of  motion  of  the  center  of  the  polar 
vortex,  have  a singularity  at  0q  = 0 , which  comes  from 
the  system  of  coordinates  used  to  describe  the  motion. 

When  the  initial  displacement  is  a polar  displacement, 
this  singularity  presents  no  problem  since  0^ (t)  ^ 0 for  all  t. 
In  the  case  of  an  initial  circle  vortex  displacement, 
the  longitudinal  component  <{>q  of  the  polar  vortex  is 
not  defined  at  t = At.  In  equations  (4.1),  the  function  f, 
corresponding  to  the  equation  defining  <J>Q  , has  a multipli- 
cative factor  of  the  form  l/(sin  0^)  . 

Numerical  experiments  were  carried  out  with  two 
different  methods  to  overcome  the  above  difficul.y: 
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(a)  At  t = 0 , displace  both  the  polar  vortex  and  a circle 
vortex,  in  such  a way  that  the  polar  vortex  displacement  is 
equal  to  a small  fi'action  of  that  of  the  circle  vortex.  The 
effect  of  the  initial  displacement  at  the  pole,  is  negligible 
compared  with  the  effect  of  the  initial  circle  vortex 
displacement  if 


where  rQ  is  the  initial  angular  displacement  at  the  pole  and 
is  the  initial  angular  displacement  of  the  circle  vortex. 
In  this  way,  Og(t)  is  nonzero,  the  singularity  is  avoided 
and  the  resulting  motion  has  the  overall  characteristics 
of  a circle  vortex  perturbation.  In  the  numerical  calcula- 
tions we  reported  in  Section  3,  values  of  | e 0/G 1 1 were  taken 
to  be  0.01,  0.02  and  in  one  case  0.1.  From  the  equations 
of  motion  given  in  Section  1,  it  is  easy  to  see  that  the 
disturbance  created  by  the  initial  displacement  of  the  polar 
vortex  is  of  the  order  of  Eg. 

(b)  If  the  velocity  vector  q of  the  polar  vortex  is 

considered  to  be  in  the  tangent  plane  (x,y)  at  the  pole, 

then  its  components  (q  ,q  ) in  that  plane  are  given  by 

x y 


(4.2) 


<*x  1 

I 

-sin  4>q  - cos  <J>Q 

' U0  ' 

q J 

cos  <Pq  - sin 

1 vo 

when  the  polar  vortex  is  considered  to  be  at  0 = 0,  <f>  = <f>Q 
and  to  have  spherical  velocity  components  (Uq,Vq). 


-90- 


On  the  other  hand,  at  t = 0,  the  polar  vortex  satisfies 
0 U'  ^0  = arbitraiy.  We  can  calculate  uQ  = a sin  , 

v0  " ‘ a00  using  the  equations  of  motion  of  Section  1. 

We  find  that 

N 

(4'3)  U0  = - 4a  sin  I Pv  ( -cos < °k+c  \ ) ) sin  (O^ciy 

K—  X 

•cos  (4>k+cuk-4)  ) 

= u cos  4)^  + v sin  <}>q 

N 

(4-4)  V0  = 7a"  sin  vtt  J.  Pv  ( _cos  ( VcBk  > > sin  (Vc(V 

K=  X 

•sin  (4>k+cnk-<j>0) 

= -u  sin  <J>  + v cos  <t>  . 

0 0 

Hence,  the  (x,y)  coordinates  of  the  center  of  the  polar 

vc  tcx  may  be  predicated  at  At  = t from 

x = At  q , 

1x 

y - 4t  qv  . 

-i 

Using  (4.2),  (4.3)  and  (4.4),  we  see  that 

qx  = -v 

qy  = u . 

If  we  map  the  coordinates  x,y  to  the  sphere  by  means 
of  a stereographic  projection, 

9 = arccos  (— 2 — ■— 2+^ 2~~ ) 

4a  + x +y 

<t>  = arctan  , 

we  find  predicted  values  for  0Q  (At)  = 0*  and  4> Q (At)  = 0*. 
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Now  by  using  a first  order  corrector  formula,  the 


polar  vortex  is  located  at  At  from 

4>0At  = 4>0°  + At  g 0i°,(J)i0) 

90At  = 90°  + At  h ( 0*,  , 0i° , 4>j_° ) , i = 1,...,N, 


where  g,  h can  be  found  using  the  equations  of  motion 
of  Section  1.  We  could  have  obtained  corrected  values 
for  0q  (At)  and  4>g(At),  accurate  to  second  order 

by  using  an  iteration  of  the  form 


u 


t+  ; 


= ufc  + 


At  f(t  + — 


At 

2 


t * 
u +u 


With  this  starting  method,  we  continued  the  original 
numerical  integration  sequence  (4.1)  until  0q  become 
smaller  than  a prescribed  value,  say  10  ‘ degrees,  at 
which  point  a similar  mapping  was  needed  to  continue 
the  integration  until  the  polar  vortex  passed  away  from 
the  pole . 

Some  integrations  were  carried  out  by  using  both  alter- 
natives (a)  and  (b) . The  results  show  no  significant 
differences.  Hence  method  (a)  was  used  primarily  owing 
to  its  simplicity. 

The  spherical  harmonic  P^ (x)  that  defines  the  stream 
function  for  geostrophic  vortices  on  a rotating  sphere, 
satisfies  the  Legendre  equation 

2 

(4.5)  (1-x2)  P"(x)  - 2 x P^(x)  - -4---  ~~  Pv(x)  = 0 . 

where 
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A 

v 


2 2 

1 f!6u'  a 

2 ( 9 h0 

- j + i A 


x 


cos  a . , 
lk 


= - cos  0k  cos  0i-  sin  0R  sin  0i  cos  ($k-<j>.)  , 

coordinates  of  the  £th  vortex  point. 

P (x)  can  be  written  [7]  : 


P (x) 


= F(j  - iX,  \ +i X , 1;  i^L) 


- 1 ^ 1 + 4X  rl-X'i  , (1+4A2) (3+4A2)  , 1-x,  2 

22  1 2 J ^2 

. (1+4A2)  ( 324-4X2)  (52  + 4A2)  ,1-^3 

i».4J.62 l~]  ••• 


= 1 + a^  + di^  + . . 


for  -1  < x < 3. 


From  the  series  expansion  of  P (x 


n+1 


= a 


n 


( (2n- 1)  2 + 4 A 2 ) 
( 2n)  2 


) , we  see  that 


a 


1 


The  squares  of  the  even  and  odd  integers  can  be 

calculated  once  and  stored.  Furthermore,  for  fixed  X, 

2 2 

we  calculated  and  stored  -2n-1)  +*A  . The  calculation 

( 2n) 

of  one  term  of  the  series  a^  requires  only  2 multiplica- 
tions and  therefore  the  recursive  calculation  of  a is 

n 

an  efficient  way  to  calculate  P (x) . Since  P (x)  has 

v v 

a logarithmic  singularity  at  x = -1 , for  values  of  the 
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argument  - 0.995  <_  x <_  -0.94,  i.e.  5°  < 20°,  about 

400  terms  of  this  series  were  necessary  to  get  a relative 

I 

accuracy  of  better  than  1/1000  for  the  value  of  P (x) . 

The  first  derivative  of  Py (x)  was  evaluated  by  using 
the  derivative  of  the  series  for  P^ (x)  while  the  second 
derivative  was  calculated  by  using  (4.5)  and  the  already 
computed  values  of  P and  P^. 

The  numerical  integrations  were  monitored  by  means  of 
the  invariant  (1.14), 


I (t) 


N 

l y cos 
k = 0 


In  general,  the  maximum  relative  change 


was  bounded 
than  5°,  E 


by 

was 


E = max 


.1-10 


-4 


t 


bounded  by 


iI(t)-KO)  | 

1 1(0)  1 

For  initial 

0.1  • 10~7 


displacements 


less 
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5 . Comparison_qf  Geostrophic  Vortices  on  the  Sphere 

Reel  t linear  Geos  trophic  Vortices  on  a Rotating  Plane . 

It  is  of  interest  to  make  a qualitative  comparison 
between  the  rectilinear  geostrophic  vortices  studied  by 
G.  K.  Morikawa  and  E.  V.  Swenson  in  [2]  and  the  geostrophic 
vortices  on  the  sphere  studied  in  this  work. 

A single  rectilinear  geostrophic  vortex  on  a rotating 
plane  is  defined  by  the  stream  function  [2] 

^0  = 2?  V K lr’r0l) 

where  y = strength  of  the  vortex 
K = 2^//gh^  , 


w = angular  velocity  of  the  plane 

'r_ro'  = disfcance  from  an  arbitrary  point  on  the  plane 
to  the  center  rQ  of  the  vortex 
Kq(x)  = Bessel  function  of  order  0. 


The  stream  function  for  a distribution  consisting  of  one 
vortex  with  strength  yQ  at  the  origin  of  the  system  of 
coordinates  and  N vortices  of  equal  strength  y,  equally 

spaced  on  a circle  of  radius  r,  centered  at  the  origin, 
is  given  by 

Y0  v N 

* ‘ 2?  V K lr-rol>  + k I K <k  |r-r.|>  [21. 

i=l  •L 

"mce  in  I2!  the  radius  of  the  circle  is  normalized  to  1, 
for  the  plane  problem  there  are  three  independent  parameters, 
(1)  N,  the  number  of  vortices;  (2)  < , scale  parameter, 

(3)  Vy,  quotient  of  strengths. 
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In  the  spherical  case,  for  a distribution  consisting 
of  one  polar  vortex  with  strength  yc  , N vortices  of  equal 
strength  m , equally  spaced  in  a circle  of  colatitude  6, 
no  boundary  condition  at  the  equator,  we  have  4 independent: 
parameters : 

(1)  N,  number  of  vortices,  (2)  6,  colatitude  of  the  circle 
vortices,  (3)  Wq/M/  quotient  of  strengths,  (4)  cp  , coriolis 
parameter. 

The  first  observation  is  that  the  spherical  problem 
depends  on  one  more  parameter  and  a direct  parametric 
equivalence  cannot  be  easily  established. 

However,  two  important  associations  between  the  stream 
functions  can  be  made: 

(1)  Peters  [3]  noted  that 


lim 

a-*°° 

where 


4 sin  vtt 
2o)  *c 


P (-cos  a.)  = - 


4tt  K0 


f~2  2~ 

K'/r  +ri" 


2r  r-^cos  ((ft-^) 


k = 


v 


*-?-?/( 


^h0 

a = radius  of  the  sphere. 


16u)2*Cp*a2  / ghQ)  - 1 


(2)  As  it  was  stated  in  Section  2,  for  small  values  of  , 
the  stream  function  for  one  vortex  in  a rotating  sphere 


tp  = 


rsin  v,  Pv (_cos  al>  < 


cos  = cos  0 cos  0^  + sin  0 sin  0^  cos  ( — 


becomes  independent  of  c^ . Furthermore, 
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1-cos 

P (-cos  ^ A log  ( 2 ) • 

On  the  other  hand,  for  < close  to  0 we  have  that 

, . I r_r0 I 

K0(  K lr-r0D  % B lo<3  — 2 » 

so  that  both  stream  functions  have  the  same  type  of 
logarithmic  behavior  for  values  of  the  argument  close 
to  their  respective  singularities. 

Since  the  two  problems  have  the  same  type  of  symmetries 
and  geometric  properties,  the  matrices  of  the  linear 
problem  have  the  same  form.  This  fact  makes  some  basic 
properties  of  the  linear  analysis  common  to  both  cases: 

(a)  zero  is  an  eigenvalue  of  multiplicity  two  and  in  some 
cases  of  multiplicity  four. 

(b)  The  equilibrium  angular  velocity  ft  is  an  eigenvalue 
of  the  linear  matrix. 

(c)  There  are  two  linear  invariants.  When  a polar  vortex 
is  present,  they  are  independent  of  the  coordinates 
of  the  polar  vortex. 

(d)  Two  distinct  families  of  eigenvalues  with  positive  real 
part  are  evident,  one  family  is  pertinent  to  the  circle 
vortex  and  the  second  family  to  the  center  vortex. 

(e)  In  particular,  for  N = 2,  uQ/u  =0  is  in  the  unstable 
range.  If  pQ/y  > 0'  the  configuration  is  exponentially 
unstable  for  all  values  of  the  parameters. 

For  N j-  3 the  configurations  are  exponentially 
unstable  for  all  values  of  the  parameters  when 

Vq/v  = C*  (Section  2) . 

-97- 


(f) 


Some  important  differences  can  be  listed  as  follows: 


Planar  Case 


[2] 


Spherical  Case 


(a)  \\Q/\i  < 0 

N = 3 is  exponentially 
stable  for  all  values 
of  < . 

4 <_  N 6 , exponentially 
stable  for  some  values 


(a)  \sQ/u  < 0 

N = 3 is  exponentially 

stable  for  some  values 

of  c and  0 . 

P 

N ^ 4 exponentially  stable 

for  some  values  of  c , 0. 

P 


Of  K 

N 7 , exponentially 

unstable  for  all  values 
of  < . 


(b)  uQ/u  = 0 

N = 2, 3, 4, 5,  exponenti- 
ally stable  for  all 
values  of  < . 

N = 6,  exponentially 
stable  for  some  values 
of  k . 


(b)  uQ/u  = 0 

N = 2,3,  exponentially  stable 

for  all  values  of  c , 0. 

P 

N > 4,  exponentially  stable 
for  some  values  of  c , 0. 


N ^ 7 , exponentially 
unstable  for  all  values 
of  < . However,  for  N=7, 
K = 0 is  a neutral 
stability  point. 


Consider  a configuration  of  three  vortices  of  equal 
strength  u equally  spaced  on  a circle  of  colatitude, 
one  polar  vortex  with  strength  , no  boundary  condition, 

and  a choice  of  Uq/u  so  that  the  circle  vortices  remain 
stationary. 

Bauer  and  Morikawa  [2aj  established,  by  using  a plane 
model,  that  the  above  system  of  vortices  can  be  used  to 
represent  the  climatological  semi -permanent  highs  whenever 

1 ’ 45  - 5'°*  These  values  of  < imply  that  the  ratio 

q,  quotient  between  the  long  period  and  short  period  of 
the  system,  takes  on  values  between  2 and  32.85. 

In  the  spherical  case,  values  of  q between  2.0  and 
10.0  are  found  at  colatitudes  55°  <_  0 < 75°,  which  is 
approximately  the  geographical  region  where  the  semi- 
permanent highs  are  observed. 

In  particular,  if  we  take  the  magnitudes  of  and  u 
in  such  a way  that  2000  Km  away  from  the  center  of  a 

circle  vortex,  the  speeds  are  about  10  [m/sec],  we  obtain: 

(a)  Planar  case 


9 = 66.5° 

< = 3.75 

q = 12 

U0/u  = - 0.081 
Short  period:  1 year 


6 = 60° 

< = 3.2 
q = 7.7 
UQ/U  = - 0.121 


short  period:  0.52  years 
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(b)  Spherical  Case 
0 = 66.5° 
c =1.0 

p 

q = 8.45 
P0/M  = -.119 
short  period: 
1.6  years 


6 = 60° 
cp  = 1-° 

q = 6.89 
Mq/P=  - .144 
short  period: 
0.4  years. 


In  Figure  5.1  we  have  plotted  |Pq/m|  = |c 
for  the  following  cases: 

(1)  Plane  mod«l,  curve  A.  The  left  end  of  the 
corresponds  to  K-  3.6,  while  the  right  end 
to  k = 2.0. 

(2)  Spherical  model  with  circle  vortices  at  0= 

The  left  end  of  the  curve  corresponds  to  c 
. . P 

the  right  end  corresponds  to  c = .2. 

(3)  Spherical  model  with  circle  vortices  at  6 = 

Left  end:  c = 1.0,  right  end:  c = .4. 

p P 

(4)  Spherical  model  with  circle  vortices  at  0 = 

Le^t  end:  c^  = 1.0,  right  end:  c = .5. 

P 


'|  versus  q 
curve 

corresponds 

70°,  curve  B. 
=1.0  while 

60°,  curve  C. 

50°,  curve  D. 
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A . Appendix 


Numerical  Results  on  the  Eigenvalue  Calculation. 


As  it  was  stated  before,  the  strength  u does  not  have 
any  effect  on  the  linear  stability  properties.  For  the 
majority  of  numerical  experiments,  it  was  chosen  to  be 


M = .15  x 1010  [— ] . 

seg 

For  this  value-,  the  speed  of  the  velocity  field  20°  away 
from  a center  of  a circle  vortex  is  between  8 [m/sec]  and 
20 [m/sec] . 

because  of  meteorological  considerations,  the  following 
choice  of  hQ,  a,  g was  made 


hQ  = 8500  [m] 

a = 6370  [km] 

g = 9.806 

seg 

Computations  were  made  for  N = 2,3,. ..,8,  c = .2, .3,..., 

P 

1.0  and  9 = 5°, 6°, . . . ,90°. 

i)  N Vortices  on  a Circle  of  Colatitude,  No  Polar  Vortex, 

No  Boundary  Condition  at  the  Equator.  (2  < N < 8) 

For  N = 2,3  we  have  exponential  stability  for  0 < 0 < 

For  N > 3,  a polar  cap  of  stability  is  characteristic  for 

each  N.  Table  A.l  shows  the  values  of  the  limit  0 of  the 

c 

stability  cap.  As  we  can  see  there,  for  fixed  N,  0 decreases 

c 
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with  increasing  cp  and  for  fixed  cp,  6c  decreases  with 
increasing  N.  Table  A. 2 shows  the  values  of  0c  with 
0^  = cos  0,  i.e.,  using  the  true  value  of  the  coriolis 
parameter. 

Taole  A. 2 Values  of  0„  tor  c = cos  0. 

c p 

Configuration  is  stable  for  0 <_  0 Q c* 


N 

6c 

4 

- 

50.4093° 

5 

37.9224° 

6 

17.5225° 

7 

9.9952° 

8 

8.4092° 

Let  w(N,e,c  ) = max  |lm  A.  I . 
p i 

This  function  will  give  us  an  overall  view  of  the  behavior 
of  the  maximum  frequency  with  respect  to  N,  0,  c . Its 

lr 

properties  are  the  following: 

1.  For  fixed  N it  is  continuous  with  respect  to  0 and 

c for  0 < 0 < tt/2  , .2  < c < 1.0. 

p - P - 

2.  lim  w(N,e,c  ) = 00  for  all  N , c . 

0^0  p p 

3.  < oj  (N , 0 , c ) for  all  N,0,  c . 

— p y 

4.  ||  (Nf0,c  ) 1 0 for  all  N,  cp  and  for  0 <_  0c- 
For  0 <_  0c  the  eigenvalues  of  L have  the  following 

distribution . 

1.  1 pair  of  0 eigenvalues; 

+ i ; 


2. 
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3.  n-2  pairs  of  imaginary  complex  conjugate  eigenvalues. 

For  0 > 0 + 0.0001°  we  have 

— c 

1.  Double  zero  root 
2 . + i SI 

3.  Eigenvalues  of  the  form  4 a^ , 4 a2 ±ib2 , a1,a2,b2  real. 

In  Figure  A-l,  che  curves  of  linear  stability  with 

respect  to  0 and  c^  have  been  drawn.  For  N fixed,  we  have 

linear  stability  to  the  left  of  each  curve. 

As  we  can  see,  when  N increases,  the  curve  becomes  more 

parallel  to  the  c axes  and  therefore  the  value  0 does  not 

P c 

depend  so  much  on  c . This  result  stems  from  the  fact 
stated  earlier:  when  0 is  small,  the  eigenvalues  are 

independent  of  cp.  A graphic  representation  of  max  | Re  Ai | 
versus  colatitude  (0)  and  Im  | A ± | , i = 1,...,N  for  different 
values  of  the  parameters  are  given. 
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= 
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i a. . N__  Vort ices  on  a Circle  of  Colatitude L No  Pol ar  Vortex  , 

Boundary  Condition  at  0 = tt/2.  (2  _<  N <_  8) 

For  N > 2,  we  have  a cap  of  stability.  In  particular 
for  N = 2,3  we  ]o  not  have  exponential  stability  for 
0<  tt/2  as  the  previous  case  but  a polar  cap.  For  N > 4, 

we  have  a polar  cap  of  stability,  but  this  cap  is  smaller  than 
that  of  (i)  for  each  N. 

In  Table  A. 3 we  have  the  limit  0 of  the  cap  of 

c 

stability  for  different  values  of  N and  c . In  Table  A. 4 we 

have  the  values  of  fl  for  c = cos  0. 

c p 

If  we  compare  the  entry  for  N 8 in  Tables  A.  2 and  A.  4 

we  can  see  that  they  are  equal  up  to  four  decimal  digits; 

th«  same  happens  in  Tables  A.l  and  A.  3 for  N = 8,  c = .9  and 

c 1.0.  This  suggests  that  when  y is  small,  the  stability 
P 


properties  of  both  cases 

are  similar, 

independent 

of 

existence  of  a boundary 

condition 

• 

We  observe  that  a . 

3 

, br  c.. 

(2. 

16) , depend 

on 

P (-cos  Y 

V 

2 

- sir.  Y 

cos 

“ki>  - 

P" ( -cos^y 

and  2 

. 2 

- sin  y 
2 

cos 

“ki>  •' 

P ' (cos  Y 

V 

- sin  Y 

cos 

“ki>  < 

P^(cos2  y 

. 2 
- sm  Y 

cos 

■Vi* 

where  the  first  pair  of  functions  comes  from  the  vortices 
on  the  northern  hemisphere  and  the  second  pair  from  the  images 
in  the  souterh  hemisphere. 


: 
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Table  A. 3.  Values  of  0 . Configuration  is  stable  for 


M 


As  y "*■  0 / the  first  pair  goes  to  + 00  respectively 
and  the  second  pair  tends  to  a constant;  hence  the  contribution 
of  the  images  is  very  small  compared  with  the  contribution 
of  the  vortices  on  the  northern  hemisphere,  and  therefore 
the  linear  stability  properties  are  the  same  for  both  cases 
when  y is  small.  From  this  we  can  say  that  for  large  N, 

the  cap  of  stability  tends  to  be  the  same  in  both  cases. 

Table  A.  4.  Values  of  0 for  c = cos  9. 

c p 

Configuration  is  stable  for  0 < 6c» 


It  is  also  worth  mentioning  that  SI  is  larger  when  we 
have  a boundary  condition;  this  is  due  to  the  increase,  in 
the  longitudinal  component  of  the  velocity,  that  comes 
from  th  - images. 
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In  Figure  A. 26,  we  have  the  curves  of  linear  stability 
with  respect  to  0 and  c . For  fixed  N,  we  have  linear 
stability  to  the  left  of  each  curve. 

If  we  compare  Figures  A.l  and  A. 26,  we  see  that  for 
N = 6,7,8,  the  curves  have  the  same  form;  if  N < 6 there 
are  noteworthy  differences.  For  N = 2 end  5,  they  intersect; 
also  note  that  for  N = 2,3,  dc  /d6  > 0 and  for  N > 3, 

dCp/dQ  < 0.  For  N =4,5,  the  curves  are  shifted  to  the 

right . 

In  general  we  can  conclude  that  the  addition  of  a 
boundary  condition  makes  the  configurations  mai?e  unstable 
in  the  sense  that  the  zones  of  stability  are  reduced. 

A graphic  representation  of  max  | Re  A^|  versus 
colatitude  and  Im  | ^ | , i = 1 , . . . ,n  , for  different 
values  of  the  parameters  are  given. 
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111  N_yorticGs  on  a Circle  of  Colatitude , 

"n-  Polar  Vortex,  No  Boundary  Condition  at  0 = j.  (2<N<8) 

Calculations  wore  made  for  N = 2,3,... ,8,  c = .2,... 1.0, 

P 

5° , 6 0 , . . . ,90  0 and  values  of  UQ/u  = c* , UQ/U  > 0, 

u0/u  ‘ 0 * Tho  general  results  can  be  summarized 
as  follows. 

(a)  ?i  = 2. 

If  lQ  M > 0 we  always  have  exponential  instability 
for  all  values  of  0,  c , un/\i . 

If  Uq/m  < 0 we  have  linear  stability  in  a belt  whose 

limits  are  0 and  tt/2. 
c 

The  table  below  gives  some  values  of  0 for 

c 

Cp  = 1-°  and  for  m0/m  = -0. 125,-0. 25,-0. 5,-1. 0,-1. 25, 

-1.5. 


!V</“ 

4 

-0.125 

i 

-0  25 

1 -0.5 

1 

0 

r— ) 

1 

-1.25 

-1.5 

.2 

87° 

83° 

74° 

46° 

11° 

nr7 

. 3 

87° 

83° 

74° 

47° 

11° 

< 5° 

.4 

87° 

83° 

74° 

48° 

12° 

< 5° 

.5 

87° 

83° 

74° 

50° 

13° 

< 5° 

.6 

87° 

83° 

74° 

51° 

18° 

< 5° 

. 7 

87° 

83° 

74° 

52° 

24° 

< 5° 

.8 

87° 

j 

83° 

75° 

52° 

29° 

| 

< 5° 

.9 

! 

87° 

83° 

75° 

53° 

32° 

l 

< 5° 

1.0 

87° 

L 

83° 

[ 

75° 

l 

54° 

L 

35° 

< 5° 

1 

Table  of  ©c  for  N = 2.  Configuration  is  stable  for  0 <0<_  tt/2. 
For  = c*  1:110  configuration  is  unstable. 
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As  we  can  see  when  |u0/u|  increases,  6c  decreases  and 
the  configuration  becomes  stable  for  q*  <•  g < 90o  0*  < 5„ 

(b)  N = 3. 

M0/u  < 0.  For  a(cp)  < pQ/p  < 0,  max  |a(c  )|  = 2.0, 

we  have  stability  in  the  whole  northern  hemisphere. 

F°r  Uq/m  i -*  a ( cp) , the  stability  region  is  reduced  to 

a cap  that  decreases  with  increasing  | m0/u | . In  the  table 

below  we  have  the  limit  6c  of  the  cap  for  different  values 

of  un/u  and  c . 

0 P 

Table  of  0c  for  N = 3.  Configuration  is  stable  for  0 < 6 < 6 


For  u0/u  - c*  we  have  stability  for  all  6,  c . 


UQ/u  > 0. 


For  0 < po/M  < b (Cp)  , max  b(cp)  = 0.75, 
we  have  stability  for  0 < 9 < n/2 . 

If  UQ/U  increases,  a belt  of  instability  appears, 
and  the  width  of  this  belt  increases  with  increasing  Vq/u. 

In  the  table  below,  we  have  the  limits  of  the 
instability  belt  01  and  02  — unstable  for  ei  < 0 < e2 
stable  otherwise , for  some  values  of  cp  and  Uq/u . 


Configuration  is  unstable  for  < o < o 
(in  degrees) 


•Sin/U 

cp^N 

1 

°i 

.0 

°2 

1 

.25 

°2 

1 

°i 

. 5 
°2 

2 

°1 

'.0 

02 

2 

61 

.5 

3 

91 

. 0 
°2 

PI 

.0 

pi 

5 

°1 

.0 

°2 

.2 

o 

0 

ii 

28 

9 

37 

9 

47 

8 

52 

8 

56 

7 

60 

7 

63 

.3 

0 

0 

10 

29 

9 

39 

9 

49 

8 

54 

8 

58 

7 

62 

7 

65 

.4 

0 

0 

10 

32 

9 

41 

8 

51 

8 

57 

8 

60 

7 

64 

7 

67 

.5 

0 

0 

10 

34 

9 

44 

8 

59 

8 

59 

8 

63 

7 

67 

7 

69 

.6 

0 

0 

10 

37 

9 

46 

8 

57 

8 

62 

8 

66 

7 

70 

7 

72 

.7 

14 

17 

10 

39 

9 

49 

8 

59 

8 

65  . 

8 

68 

7 

72 

7 

74 

.8 

13 

22 

10 

42 

9 

51 

8 

62 

8 

67 

8 

71 

7 

75 

7 

77 

. .9 

12 

25 

10 

43 

9 

53 

8 

64 

8 

69 

7 

73  i 

7 

77 

7 

79 

1.0 

12 

28 

10 

44 

1 

9 

54 

8 

1 

65 

1 

8 

J 

71 

l 

7 

75 

7 

79 

7 

81 

(c)  N > 4 . 

Uq/m  < 0,  for  small  values  of  |pg/y|  we  have  a polar 
cap  of  stability  which  decreases  with  increasing  values 
°f  | Uq/u | • For  Pq/m  < -1.0,  the  configuration  is  unstable 
in  the  range  5°  < 6 < 90°.  Given  below  are  the 

values  of  the  limit  6c  of  the  cap  for  N = 4, 5, 6, 7. 


\ 


For  Mq/u  = c*  the  configuration  is  exponentially 
unstable. 
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For  ^q/v  = c*  we  have  exponential  instability 
for  N = 4, 5 , € , 7 , 8. 
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If  Mq/p  > 0,  we  have  a cap  of  stability. 

As  |Jq/|j  increases,  this  cap  becomes  larger,  reaching 
a point  where  vse  have  stability  for  0 <_  6 <_  tt/2.  Then  a 
belt  of  instability  appears.  With  increasing  Pq/P/  the 
belt  of  instability  increases.  In  the  tables  given  below, 
we  have  the  limit  0C  of  the  cap  and  the  limits  6^,  ©2  of 
the  belt  for  some  values  of  N,  p0/p , c . 

\J  P 
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N = 5 


Values  of  0 . 

c 


Configuration  is  stable  for  0 < 9 < 0 

— — c * 


— 

.25 

1 -5 

i 

| 1.0 

1 

1.5 

1 

3.0 

.2 

| 

53° 

1 

61° 

| 

73°  ; 

82° 

90° 

•3 

53° 

61° 

73° 

83° 

90°  i 

.4 

53° 

1 

61° 

74° 

83° 

90° 

•5 

53° 

62° 

75° 

84° 

90° 

.6 

1 

52° 

62° 

76° 

84° 

90° 

•7 

1 

52° 

63° 

76° 

84° 

90° 

•8 

52° 

64° 

77° 

85° 

90° 

.9 

53° 

64° 

78° 

85° 

90° 

1.0 

J 1 

53° 

l 

65° 

l 

78° 

L 

85° 

L 

90° 

1 

Values  of  0^  02. 

Configuration  is  unstable  for  01<0<02;  stable  otherwise. 


XV1 
!cp  X 

4 

.5 

5 

.0 

10 

.0 

11.0 

9i 

92 

ei 

92 

6! 

02 

! 9? 

J 

13° 

j 15° 

11° 

21° 

9° 

38° 

n — 

i 8° 

i 

40° 

. 3 

i 

j _ 

11° 

20° 

9° 

37° 

8° 

38° 

.4 

1 _ 

11° 

18° 

go 

350 

8° 

37° 

.5 

- 

12° 

17° 

go 

34° 

8° 

35° 

.6  1 

( 12° 

16° 

go 

32° 

j 8° 

33° 

.7 

- 

1 1 2 0 

14° 

go 

30° 

8° 

31° 

.8 

- 

- 

13° 

13° 

9° 

29°  i 

• 

8° 

30° 

.9 

- 

- 

- 

go 

27° 

| 

8° 

28° 

1.0 

L 

LLi 

' | 

LJ 

I 

L 

8° 

26° 

1 

L 

i 

8° 

; 

27“  ! 

L 
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values  of  0 

c 

Configuration  is  stable  for  0 < 0 < 6 

— — c 


►— * 
• 

o 

1.5 

3.0 

5.0 

55° 

64° 

1 

81° 

90° 

55° 

64° 

81° 

90° 

55° 

64° 

81° 

90° 

55° 

64° 

80° 

90° 

54° 

64° 

80° 

90° 

54° 

64° 

79° 

90° 

54° 

64° 

79° 

88° 

54° 

63° 

78° 

87° 

54° 

63° 

J 

77° 

I 

86° 

Values  of  0^ , 02 


Configuration  is  unstable  for  O^Q^;  stable  otherwise 


M 


7.5 

8 " 

61  e2 

6l|  e2 

12°  17° 

11°  20° 

12°  16° 

11°  19° 

O 

• — 1 

o 
CM 
• — < 

11°  17° 

_ 

11°  16° 

- 

11° | 14° 

__ 

12°  12° 

- 

1 - 1 - 

10°  23° 
10°  22° 
10°  20° 
10°  18° 
10°  17° 
10°  15° 


11.0 

ei|  e2_ 

9°  29° 
9°  27° 
9°  25° 
9°  23° 
9°  22° 
9°  20° 
9°  18° 
9°  17° 
9°  16° 


N 


7 


Values  of  0 

c 

Configuration  is  stable  for  0 < 0 < 6 . 

- — c 


A Graphic  Representation  of  max  I Re  X^  | and  Im  | X^  | for 
Different  Values  of  the  Parameters  are  Given 


Figure 

Number 


Con  tents 


Plot  of 


max 

| Re 

xii 

0 

3 

• 2 f • • • 9 

1.0 

max 

| Re 

Ail 

9 

3 

• 2 f 0 • • r 

1.0 

max 

| Re 

Ail 

e 

3 

• 2 f ...  t 

1.0 

max 

1 Re 

x.  1 

0 

3 

• 2 t • • • t 

1.0 

1 

Im 

Pi 

e 

3 

1.0 

max 

| Re 

Ail 

e 

4 

• 2 1 0 % • t 

1.0 

max 

| Re 

Ail 

e 

4 

•2 f 0 0 0 t 

1.0 

max 

| Re 

Ail 

e 

5 

• 2 ! • • • t 

1.0 

max 

| Re 

Ail 

e 

5 

0 2 1 • • • t 

1.0 

max 

| Re 

Ail 

e 

5 

• 2 t . • • t 

1.0 

max 

| Re 

X . I 

e 

6 

• 2 ! • • • t 

1.0 

1 1 

Number 


Contents 


For 
N = 


Numoer 

Plot  of 

Versus 

N 

A.  38 

max  |Re  A | 

0 

6 

A.  39 

max  1 Re  A . 1 

l 1 

0 

6 

A.  40 

max  1 Re  A . 1 

l ' 

0 

6 

A.  41 

max  | Re  A . 

l 

0 

7 

A.  42 

max  1 Re  A . 1 

' l ' 

0 

7 

A. 43 

max  1 Re  A . 1 

1 l 1 

0 

7 

A. 44 

max  |Re  A . ( 

1 x 1 

0 

8 

A.  45 

max  | Re  A . 1 

l 1 

0 

8 

A.  46 

max  1 Re  A . 1 

l 

0 

8 

A. 47 

max  | Re  A . 1 

l 

0 

8 

au»il  r-xe 


IC.atllUOf 


Figure 

lh« 

HUO/KJs  - SOOO 


.31 


Figure  A. 32 


* n=s 


4 1 HUO/HUs  .sooo 


OKI  iWi  XU 


COlfii  ITJO! 


CCH  AT  | T UOC. 


Fiaure  A. 44 


1.5QO 


njo/Wz  i.ooro 


t*  miruot 


5 


N Vortices  on  a Circle  of  Colatitude,  One  Polar  Vortex. 
Boundary  Condition  at  0 = tt/2.  (2  < N < 8) 

(a)  N = 2.  Uq/p  > 0. 

We  have  a belt  of  stability,  the  limits  and  ©2  of  the 
belt  are  given  below  for  some  values  of  and  UQ/u- 


Configuration  is  stable  for  <_  0 £ e 
(in  degrees) 


c \ 
p \ 

0_^ 

125 

• 

25 

5 

1 

.0 

1 

.5 

3 

.0 

— 

4 

. 0 

1 

°1 

62 

61 

62 

61 

02 

91 

°2 

°1 

l02 

-! 

°1 

j02 

61 

-1 

! 

.2 

65 

66 

67 

70 

70 

76 

72 

82 

72 

84 

71 

— 

>85 

63 

>85 

rj 

.3 

65 

66 

67 

71 

69 

76 

71 

82 

72 

84 

71 

>85 

63 

>85 

! 

.4 

65 

66 

67 

71 

69 

70 

71 

82 

72 

84 

72 

>85 

70 

>85 

i 

: 

.5 

65 

67 

67 

71 

69 

77 

71 

82 

72 

84 

72 

>85 

71 

>85 

j 

.6 

65 

67 

67 

72 

69 

77 

70 

81 

71 

84 

72 

>85 

71 

>85 

j 

.7 

65 

68 

66 

72 

68 

77 

70 

81 

71 

83 

72 

>85 

71 

>85 

| 

.8 

65 

69 

66 

73 

68 

77 

70 

81 

70 

83 

71 

>85 

71 

>85 

1 

.9 

64 

69 

66 

73 

68 

77 

69 

| 

81 

70 

83 

71 

>85 

7lj 

>85 

1.0 



J 

64 

70 

J 

66 

__l 

74 

67 

78 

_i 

' 69 

1_J 

81 

L 

1 70 

J 

83 

L 

71! 

l 

>85 

l 

71 

| | 

>85  1 


j 

For  3.5  <_  Uq/ u < 4.0  another  belt  appears  between 
60°  and  70°  colatitude.  The  values  for  this  belt  are  not 
tabulated.  It  can  be  seen  in  Figure  A. 69  a. 
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Wg/W  < 0* 

For  -1.0<  yo/M  < 0 «e  have  instability.  For  uQ/u  = -1-25 

, . n /n  < -1.5  we  have 

L belt  of  stability  appears  and  for  -2  _ mq/  _ 

i cap  of  stability. 

The  limits  and  92  of  the  belt  of  stability  and 
the  limit  9C  of  the  cap  are  given  below  for  various  values 

Of  Cp  and  y q/m  - 

• cfnhlP  for  0 < 6 < 09?  and  e*i  9 - ec  ' 0 <5° 

Configuration  is  stable  fo  i - - 2 


XJo7*7- 

C \ 

D \ 

-1. 

°1 

25 

6 2 

-1.5 

6c 

-2.0 

e 

c 

...  k— 

.2 

11° 

47° 

57° 

65° 

.3 

11° 

45° 

55° 

63° 

p 

12° 

43° 

53° 

61° 

.5 

13° 

42° 

O 

H 

m 

5R° 

.6 

18° 

40° 

50° 

57° 

.7 

26° 

49° 

55° 

.8 

32° 

37° 

48° 

54° 

.9 

unstable 

48° 

54° 

1.0 

unstable 

-L 

48° 

53° 

1 

For  U0/U  < -2.5  , this  cap  is  broken  into  a smaller 
cap  and  a belt  of  stability.  In  the  tables  below  we  have 
the  limits  of  the  cap  and  belt  for  different  values  of  cp 

and  Mg/M  • 


-142- 


Uq/m  = -3.5 


Configuration  is  stable  for  0.^  <6  <_  e 


VQ/\i  = -4.0 

Configuration  is  stable  for  9X  I 9 i ®2' 
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- — — — 


■tMt 


Pq/m  = -4.5 

Configuration  is  stable  for  6^.  <6  <_  e 


N 


3 


(a)  yQ/y  < 0. 

In  this  case  we  have  a cap  of  stability. 


Given  below  are  some  values  of  0 . 

c 


Value  of  0 . Configuration  is  stable  for  0 < 0 < 0 . 

c J - - c 


(b)  \,Q/\i  > 0. 


For  0<  Vq/u  < .85,  we  have  a cap  of  stability  defined 

by  0 < 0 . In  the  table  below  we  have  the  values  of  0 
— c c 

for  Uq/u  = .5,  .75,  .85. 

Values  of  0c«  Configuration  is  stable  for  0 £ 0 <_  0 . 


v\p0/u 

V 

_cp  \ 

0.25 

.5 

.75 

. 85 

.2 

46° 

47° 

48° 

49° 

.3 

46° 

47° 

48° 

48° 

. 4 

46° 

47° 

48° 

48° 

.5 

46° 

47° 

48° 

48° 

.6 

46° 

47° 

48° 

48° 

.7 

46° 

47° 

47° 

48° 

.8 

46° 

47° 

47° 

47° 

.9 

46° 

47° 

47° 

47° 

1.0 

46° 

47° 

— — . 

47° 

47° 

For  larger  values  of  yQ/u,  .85  < uQ/u  < 1.5,  this 
cap  of  stability  is  broken  into  a smaller  cap  and  a belt 
of  stability. 

In  the  table  below  we  have  the  limits  of  the  smaller 
cap  and  the  belt  of  stability  for  values  of  Pq/u  = 1.0  and  1.5. 


Configuration  is  stable  for  0 < 0 < 0 and  0 < 0 < fi 

- - c 1 - - 2 ’ 


CP 

i~ 

= 1 

0 

J0/M  = 1 

.5 

[cap  lim 

0 

c 

limi  ts 
°1 

of  belt 
’ °2 

cap  lim 

0 

c 

limits 
°i  - 

of  belt 
6 2 

.2 

49  0 

- 

- 

— 

8° 

43° 

52° 

*3 

49° 

- 

- 

8° 

45° 

51° 

' .4 

49° 

- 

- 

8° 

47° 

50° 

■5 

49° 

- 

- 

8° 

- 

- 

.6 

48° 

- 

- 

8° 

- 

.7 

13° 

O 

00 

rH 

48° 

8° 

- 

— 

.8 

, 

12° 

23° 

O 

00 

■^r 

8° 

- 

_ 

.9  | 

11° 

27° 

47° 

8° 

- 

_ 

_ °-  i 

1 

O 

rH  L 
rH 

29° 

47°  | 

8° 

- 

: i 

The  belt  disappears  for  increasing  values  of  u /v 
and  for  uQ/\i  = 2.0,  2.5;  we  have  again  a cap  of  stability 


Values  of  0 . Configuration  is  stable  for  0 < 0 < 0 


I 


Increasing  Mq/m,  the  new  cap,  is  broken  into  a smaller 
cap  and  a belt  of  stability.  Given  below  are  some  values  of 
the  limit  0 of  the  cap  and  the  belt  for 

different  values  of 

Configuration  is  stable  for  0 < 0 < 0c  ; i 6 1 02* 


N > 4 


PQ/p  < o. 

A cap  of  stability  is  again  common  for  n _>  4 . The  limits 

of  the  cap,  0c  , are  given  below  for  N = 4,5,6, 7,8  and 

ct  = *2,  .3,  ...,  1.0  for  different  values  of  p /ij  . 

F 0 

N = 4.  Values  of  0 . Configuration  is  stable  for  0 < 0 < 0 

c — — c. 


c NV 

- .i 

-0.25 

— 

- 0.5 

O 

H 

1 

.2 

39° 

33° 

10° 

<5° 

.3 

39° 

32° 

90 

< 5°  , 

.4 

38° 

31° 

9° 

< 5° 

.5 

38° 

29° 

90 

< 5° 

.6 

37° 

28° 

8° 

< 5° 

•7 

36° 

26° 

8° 

< 5° 

. 8 

35° 

24° 

8° 

< 5° 

.9 

34  c 

22° 

8° 

< 5° 

L1:0 

32° 

20- 

70 

< 50  1 

M = 5 

Values  of  0 . Configuration  is  stable  for  0 < 0 < 0 . 


M0/M  > o. 


Here  we  have  a cap  of  linear  stability,  as  \iQ/\ j 

increases,  this  cap  is  broken  into  a smaller  cap  and  a belt 

of  stability.  Below  are  given  the  limits  9c  of  the  cap, 

the  limit  01  of  the  smaller  cap  and  the  limits  02,©3 

of  the  belt  for  different  values  of  N,  un/u  and  c 

0 p* 

H - 4.  Values  of  6^  Configuration  is  stable  for  0 < 0 < e 


.25 

.5 

T' 

1.0 

pm 

p \ 

.2 

47° 

50° 

52° 

54° 

.3 

47° 

50° 

52° 

54° 

.4 

47° 

50° 

52° 

54° 

.5 

48° 

50° 

52° 

54° 

.6 

48° 

50° 

52° 

54° 

.7 

48° 

50° 

52° 

54° 

.8 

49° 

51° 

52° 

54° 

.9 

49° 

51° 

52° 

54° 

1.0 

50°  | 

51° 

52° 

54° 

Configuration  is  stable  for  0 < 0 < 0^  02<0£03;  unstable 

otherwise 


3.0 


5.0 


Configuration  is  stable  for  0 < 0 < Qy  02  < o < 0 

unstable  otherwise. 


J~ uQ/u  | 

• - 

J j 

capi  belt  1 

E- 

j£s£ 

Lh. 

5.  0 
bolt 

bl  ,J3 

!J— 

10 . 0 1 

11.0 

1 

c \ 
1 — ^ 

0 

0 

- — 4 

1 e3  I 

' i cap 

II  e, 

belt  Jj 
• 0 2 1 3 ' 

[cap 

°1 

i belt  _ 
2!  3 

T 

.2 

I 1 

12- 

16° 

63° 

|io° 

2 i o 

'64° 

8° 

39°  69° ■ 

7 J 

39° I 70° 

1 

• 3 

62° 

• 

— 

10° 

1 20° 

63° 

8° 

3 7°  b9° 

7° 

. 1 

38°  70° 

•4 

62° 

• 

— 

1 10°l 

19° 

63° 

8° 

36*  69° 

n o 
/ 

37°  6^ 

• 5 

62° 

• 

• 

1 1 0 1 

18° 

63  J 

8° 

34°  68° 

7° 

35°  69° 

. 6 

62° 

• 

111° 

17° 

62  3 

8° 

33°  68*  , 

7° 

34  0 ) 68° 

.7 

61° 

11° 

15° 

62° 

8C 

31°  67° 

7° 

32°  1 68° 

• 8 

61- 

•• 

— 

12° 

14° 

62° 

8° 

29° ! 67e 

7° 

31°  67° 

I 

.9 

61° 

— 

- 

|62°  1 

~ | 

- 

8- 

28°  66°  , 

7° 

29°  67° 

1 1 • 0 1 i 

61°  , 
| 

- | 

. l 

_n 

Liili 

-1J 

1 7°i 

27°  J 66“ j 
1 — I LL 

7°, 

2 8 0 | 66  0 | 
— L J 

' « u 1 L L 

N = 6.  Values  of  0C.  Configuration  is  stable  for  0 < ! < 8 


fx  \u0/u 

1 c 

H — 

.25 

. 5 

1.0 

1.5 

3.0 

5.0 

.2 

1 

340 

40° 

470 

52° 

59° 

63° 

. 3 

33° 

40° 

48° 

52° 

59° 

63° 

.4 

33° 

39° 

48° 

53° 

59° 

63° 

.5 

32° 

39° 

48° 

53° 

59° 

63° 

.6 

31° 

390 

48° 

53° 

59° 

62° 

.7 

31° 

390 

48° 

53° 

59° 

62° 

.8 

30° 

38° 

490 

53° 

59° 

62° 

.9 

29° 

38° 

49° 

54° 

59° 

62° 

1.0 

1 

28* 

38°  490 

L__ L 

54° 

l 

59° 

62° 

I 
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Configuration  is  Stable  for  0 <_  0 < 0^;  0^  < 0 < 

Unstable  otherwise. 


7.5 

cap  belt 

61  ^2 1 93 

11°  18°  67° 

11°  17°  66° 

11°  15°  66° 


10.0 

cap | belt 

61  V1  e3 
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7.  Values  of  0 Configuration  is  stable  for  0 < 0 < 
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Configuration  is  Stable  for  0 <_  0 <_  0 , 


°2  - 6 1 °3  ' 


\u0/u 


Unstable  otherwise. 
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N = 8.  Values  of  0 . Configuration  is  stable  for  0 < 0 < 0 , 
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Graphic  Representation  of  max  | Re  and  Im  | | versus 
is  given  for  different  values  of  the  parameters. 
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WUO/lHi 


n>c/nf=  .5000 


WUO/Wr  1.0000 


By  using  the  tables,  we  can  construct  the  curves  of 

linear  stability  which  separate  the  unstable  from  the 

stable  regions,  in  the  plane  0)  for  each  N and  c . 

In  Figures  A.69a,b,c  and  A. 70,  the  curves  are  drawn  for 

N = 2 , 8 , c = .2  with  and  without  boundary  conditions 
P 

respectively . 

Comparing  the  figures,  we  can  see  that  the  effect  of 
the  boundary  condition  is  to  reduce  the  zones  of  stability. 
Another  feature  is  that  for  smaller  colatitudes,  N >_  4 , 

6 <_  40°,  Mq/u  > 3,  the  zones  of  instability  are  almost 
the  same  with  and  without  boundary  conditions,  which  means 
that  in  the  region  mentioned,  the  stability  properties  do 
not  depend  on  the  existence  of  the  boundary  conditions. 

Also  for  small  colatitude,  the  coriolis  parameter  does  not 
play  a significant  role  in  the  stability  properties.  Both 
phenomena  were  present  when  considering  the  cases  v.’ith 
no  polar  vortex. 

Note  that  for  N = 2,  with  and  without  boundary  conditions, 
U q/u  = 0 is  in  the  zone  of  instability. 

For  the  cases  with  polar  vortex,  two  families  of  eigen- 
values are  evident,  one  family  is  pertinent  to  the  circle 
vortices  and  yield  the  left  branches  of  the  stability 
curves,  the  other  family  is  pertinent  to  the  polar  vortex 
and  yields  the  right  branches  of  the  stability  curves. 

For  the  case  of  no  pole,  and  no  boundary  condition,  we 
had  ft  = 0 at  6 = tt/2  , and  this  introduced  another  pair  of 
zero  eigenvalues. 
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When  a pole  is  added,  ft  f 0 at  6 = rr/2,  and  the 
singularity  is  removed  from  9 = tt/2  by  the  effect  of  the 

polar  vortex  and  placed  at  the  value  of  6 such  that 

* 

Uq/U  = c (0,N,Cp) • 

When  there  is  no  polar  vortex,  the  effect  of  the 
coriolis  parameter  can  be  summarized  as  follows: 

(a)  No  boundary  condition: 

For  a given  configuration,  increasing  the  coriolis 
parameter  cp  , the  polar  cap  of  stability  decreases 
for  all  N. 

(b)  Boundary  Condition 

For  a given  configuration,  increasing  the  coriolis 
parameter,  the  polar  cap  of  stability  increases  for  N = 2,3 
and  decreases  for  N > 3. 

When  a pole  is  present  there  is  no  such  simple  rule. 

For  example,  for  N = 2,  no  boundary  condition, 

Pq/m  < 0 , we  have  a belt  of  stability,  if  we  increase  c , 

P 

the  belt  decreases.  For  the  same  case  but  yQ/y  > 0,  the 
belt  of  instability  increases,  i.e.  stability  decreases, 
when  c increases.  For  N - 4,  no  boundary  condition, 

It 

Mq/u  < 0,  the  polar  cap  of  stability  decreases  with  increasing 

cp.  For  N = 5,  boundary  condition,  Vq/v  - .25,  the  cap  of 

stability  increases  with  increasing  c but  for  N = 6, 

P 

boundary  condition,  yQ/y  = .25,  the  cap  of  stability  decreases 


Figure  A. 70.  Curves  of  Linear  Stability,  N = 2, 3 , 4 , 5 , 6 , 7 , 8;  = .2. 

Polar  vortex,  no  boundary  condition.  S:  exponentially  stable  side; 
U;  exponentially  unstable  side.  No  calculations  for  9 < 5°. 


2. 


Numerical  Solutions  of  the  Nonlinear  Equations  of 
MgJr-i011  for  Confix  orations  with  No  Polar  Vortex, 

Both  with  and  Witl  ut  Boundary  Condition. 

In  Section  3,  we  havj  discussed  configurations  with 
polar  vortex  and  given  a description  of  the  motions  in  the 
zones  of  exponential  stability  and  instability. 

For  configurations  without  polar  vortex,  we  have  seen 
that  if  cx  = 0,  there  are  no  linearly  increasing  solutions; 
whereas  C1  / 0 implies  that  the  longitudinal  components 
increase  linearly  with  time. 

Configurations  consisting  of  four  vortices  on  a circle 
of  colatitude  are  studied  in  this  section.  N = 4 was 
chosen  since  it  is  the  smallest  number  of  circle  vortices 
with  configurations  in  the  exponentially  stable  and  unstable 
zones  for  both  cases,  with  and  without  boundary  conditions. 

In  the  figures  to  be  given  later  on,  we  have  not  plotted 
4^(0,  k = 1,...,4  but  instead 

Dl(t)  = % - (<*»!  (t)  + (N-l ) ^ ) 

Dk(t)  = $k(t)  - Uk.x(t)  + ~)  , k = 2,3,4. 

Dk(t)  gives  a measure  of  the  relative  oscillations, 
in  the  <j>  direction,  between  the  (k-l)st  and  k-th  vortex. 
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(i)  4 Voxtices  with  Equal  Strength  \i , Equally  Spaced  on 

a Circle  of  Colatitude  6 , No  Boundary  Condition  at 
the  Equator. 

Here  we  study  three  cases  for  the  coriolis  parameter 

c = • 7 and  one  case  for  the  coriolis  parameter  c = .5. 

P P 

The  linear  analysis  establishes  that  zero  is  an  eigenvalue 
of  multiplicity  two  for  all  values  of  the  parameters. 

In  particular  for  c = .7,  we  find  that  for 

ir 

0 <_  49.6864°,  the  configuration  is  exponentially  stable; 
whereas  it  is  exponentially  unstable  for  0 :>  49.6865°. 

For  c = .5,  the  configuration  studied  is  in  the 

IT 

exponentially  unstable  region. 

(a)  0 = 45° , cp  = • The  configuration  is  exponentially 

stab  le . 

The  experiments  were  carried  out  with  initial  displace- 
ments ef^,  ea^  between  +1°  and  + 30°.  The  results  show 
that  all  the  resulting  motions  are  linear  stable  motions, 
even  though  in  some  ca^s,  the  amplitudes  of  the  resulting 
oscillations  are  bounded  by  3/2  of  the  initial  disturbance. 
This  is  particularly  true  when  the  initial  displacement 
is  in  the  longitudinal  direction. 

For  initial  displacements  in  the  0 direction,  the 
circle  vortices  tend  to  oscillate  about  a new  position  north 
or  south  of  the  equilibrium  circle  of  colatitude  depending 
on  the  sign  of  e8„.  This  phenomenon  has  the  effect  of 
changing  the  angular  velocity  of  rotation. 
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ths  table  given  below,  we  have  listed  nonlinear 
periods  of  rotation,  for  different  types and  magnitudes  of 
the  initial  displacement. 


Configurations  at  45°  Colatitude,  c = .7,  p = .15*10 

P 

Equilibrium  Period:  677.66  [hours] 


Nonlinear  Period  Initial  Displacement 

(average) 


707.344  [hours] 

E31 

= 

1° 

722 . 894 

£31 

= 

2° 

737.723 

= 

3° 

770.285 

e3l 

= 

5° 

648.181 

= 

-3° 

636 .607 

e3i 

= 

-5.0 

677.607 

Ea, 



1° 

1 

669.103 

e a1 

= 

5° 

682.614 

ea. 

_ 

-1° 

1 

686.335 

ea. 

= 

-5° 

1 

The  change  in  the  periods  is  clearly  strong*  r for  a 
colati tudinal  initial  displacement  / 0,  see  Section  3.1) 
The  table  given  includes  nonlinear  effects,  but  for  small 


initial  displacement,  we  should  expect  that  most  of  the 
change  can  be  accounted  for  by  the  linear  theory. 

Figures  A. 71,  A. 72,  A. 73  show  some  examples  of  motions  for 


e S-j_  = 


eBi  = 


= -1°,  respectively, 


1% 


-2°, 


4 


(b)  9 - 49°,  cp  - .7.  The  configuration  is  0.6864' 


away  from  the  critical  value,  inside  the  linear  stable  zone. 

For  colati tudinal  initial  displacements,  the  motions 
present  the  same  characteriz tics  as  the  previous  configuration. 

When  the  initial  disturbance  is  in  the  longitudinal  direc- 
tion, the  situation  changes  radically.  The  amplitudes  of 
the  oscillations  of  the  colatudinal  components  are  about 
4 times  larger  than  the  amplitude  of  the  initial  displace- 
ment, while  D^(t),  k = 1,...,4,  remains  bounded  by 
the  displacement. 

If  we  decrease  the  initial  displacement,  the  amplitudes 
decrease  in  approximately  the  same  amount,  and  hence  the 
motion  can  be  considered  as  a stable  linear  motion.  In 
this  region,  nonlinear  effects  start  to  enter  into  the 
motion.  Some  runs  were  made  for  6 = 48°.  The  results  show 
that  the  nonlinear  effects  tend  to  disappear. 

In  Figure  A. 74  we  have  one  configuration  under  the 


longitudinal  perturbation  ec^  = 0.1 


The  change  in  the  angular  velocity,  presents  similar 


characteristics  as  in  the  previous  case  whenever  c^  ^ 0 


or  c^  = 0 (see  3.1). 


Configuration  at  49°,  c = .7,  y = .15*lu 

P 


10 


Equilibrium  period  948.983  hours 
Nonlinear  Period  (Average)  Initial  Displacement. 
949 . 544  hours 
948. 918 
1022.225 
1051.566 


ea^  = 0.1° 


Eax  = - 0.5° 


e3x  = 0.5° 


= - 0.5° 
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(c)  i:  8 -50.41  °,  c^  = •?•  The  configuration  is  = 0.7° 


away  from  the  critical  value,  inside  the  exponentially 
unstable  zone. 

Figure  A. 75  shows  a nonlinear  periodic  motion  with 
departure  away  from  the  equilibrium.  During  the  4000  hour 
integration,  the  motion  is  still  a nonlinear  periodic 
motion  but  eventually  it  will  depart  away  from  the  neighbor- 
hood of  the  equilibrium  position. 


ii:  8 - 60°,  c - .5.  The  configuration  is  inside  the 

IT 


exponentially  unstable  zone. 

Experiments  were  carried  out  only  with  initial 


displacements  e31#  cc^  between  -2.0°  and  2.0°. 


All  the  motions  obtained  were  nonlinear  periodic 
motions  with  departure  away  from  the  equilibrium  position. 
Figures  A.  76,  A.7"1  show  the  resulting  motions,  under  a 
positive  = 0.125°  and  negative  e3^  = -0.25°  colatitudi- 
nal  initial  displacement,  in  the  nonlinear  periodic  stage. 

From  the  overall  nonlinear  results  obtained,  we  can 
conclude  that  the  linear  analysis  gives  a good  approximation 
to  the  motion . Furthermore,  the  zones  of  exponential 


stability  and  instability  on  the  (0,c  ) plane  are  well 

P 


defined  by  the  linear  theory  for  N = 4,  in  the  sense  that 
whenever  the  linear  theory  establishes  an  exponentially 
stable  or  unstable  configuration,  it  is  confirmed  by  the 
nonlinear  calculations;  except  in  a neighborhood  of  the 
boundary  of  the  linear  stable  zone.  This  neighborhood, 


for  N = 4,  c = .7,  is  about  2°. 
P 
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List  of  Figures 

Each  of  the  figures  (A. 71)  - (A. 77)  consists  of  5 
subfigures  labeled  a,b,c,d,e.  The  subfigures  are  plots  of 
the  following  quantities: 

a:  ei(t) 
b:  02(t) 

c:  e3(t) 

d:  04(t) 

e:  Dk(t)  , k = 1,2 ,3, 4. 


For  example,  Figure  A. 72  c represents  the  component  e (t) 
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Figure  A. 75  e 
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(ii)  4 Vortices  with  Equal  Strength  m,  Equally  Spaced  on  a 
Circle  of  Colatitude  6,  Boundary  Condition  at  the 
Equator . 

Here  we  study  three  cases  for  the  coriolis  parameter 

0^  = .8.  According  to  the  linear  analysis  of  section  2, 

zero  is  an  eigenvalue  of  multiplicity  two.  Furthermore, 

for  0 <_  41.7070°,  the  configuration  is  exponentially  stable 

whereas  for  0 >_  41.7071°,  it  is  exponentially  unstable. 

Numerical  integrations  were  carried  out  for  values  of  0 

in  the  exponentially  stable  and  unstable  zone. 

(a)  0 = 45°,  c = .8.  The  configuration  is  exponentially 

P 

stable . 

Numerical  integrations  were  conducted  for  initial 
displacements  and  ea^  between  -20.0°  ana  20.0°.  The 

results  show  that  all  the  motions  are  linear  stable  motions. 
For  all  the  experiments  D^(t),  k = 1,2, 3, 4,  is  bounded 
by  3/2  of  the  initial  amplitude  while  0^{t),  k = 1,2, 3, 4, 
by  the  initial  amplitude. 

When  the  initial  displacement  is  in  the  0 direction, 
the  circle  vortices  tend  to  oscillate  about  a new  circle  of 
colat !.tude , north  or  south  of  the  equilibrium  position, 
depending  on  the  sign  of  e8^. 

This  behavior,  as  in  the  case  with  no  boundary 
condition,  effects  the  angular  velocity  of  rotation. 

In  the  table  below  we  have  listed  the  nonlinear  periods 
of  rotation  for  different  initial  displacements. 
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Configuration  at  35° , c - .8,  y = .15*10 
Equilibrium  period  327.573  hours 

(Average)  Initial  Displacement 


331.695  hours 

eB. 

X 

X 

.5 

456.822 

eB  ^ 

= 

15° 

323.508 

e&1 

= 

-0.5° 

221.274 

eBi 

= 

-15° 

327.525 

eal 

= 

.5° 

308.484 

eal 

= 

15° 

327.480 

ea^ 

= 

-1.0° 

308.153 

ea^ 

= 

-15.° 

The  change  in  the  periods  of  rotation  is  clearly  stronger 
for  a colatitudinal  displacement  (c^  t 0,  Section3.1).  The 
table  given  includes  nonlinear  effects,  but  for  small  pertur- 
bation, most  of  the  change  can  be  accounted  for  by  the  linear 
theory. 

In  Figures  A. 78,  A. 79,  A. 80  we  have  some  results  on 
the  nonlinear  calculations  for  eB^  = 0.5°,  eB^  = -15°, 
ect-^  = -15°  respectively. 

(b)  9 = 41°,  cp  = .8.  The  configuration  is  0.7170°  away 

from  the  critical  value,  inside  the  exponentially  stable  zone. 

For  a colatitudinal  displacement,  the  motions  present 
the  same  characteristics  as  the  previous  configuration. 

For  a longitudinal  initial  displacement,  the  amplitudes 
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of  Oj^t),  k = 1,2, 3, 4,  are  about  four  times  the  initial  pertur- 
bation, but  as  ca^  decreases  in  magnitude,  the  amplitudes 
decrease  in  the  same  amount.  In  this  region,  nonlinear 
effects  enter  into  the  motion. 

O 

Some  experiments  were  conducted  for  6 = 40  . They  show 

that  the  nonlinear  effects  that  produce  large  amplitudes  in 
the  colati tudinal  components  tend  to  disappear. 

(c)  0 = 50°,  c = . 8 . The  configuration  is  in  the 

exponentially  stable  range. 

The  experiments  were  done  using  perturbations  in  the 
positive  and  negative  colatitudinal  and  longitudinal  directions 
with  a maximum  amplitude  of  0.5°. 

For  all  the  calculations,  the  motions  are  nonlinear 
periodic  motions  in  the  deparature  away  from  the  equilibrium 
position.  Figures  A. 81,  A. 82  and  A. 83  present  some  motions 
in  the  exponentially  unstable  range  for  = 0.25°, 
ea^  = 0.25°,  eB^  = 0.1°  respectively. 

From  the  results  of  the  numerical  integration  of  the 
nonlinear  equations  of  vortex  motion,  we  can  state  that  the 
linear  analysis  gives  a good  approximation  to  the  motion. 
Furthermore,  the  zones  of  exponential  stability  and 
instability  on  the  (0,c  ) plane  are  well  defined  by  the 
linear  theory  for  N = 4,  except  in  a neighborhood  of  the 
boundary  of  the  exponentially  stable  range.  This  neighbor- 
hood is  about  1.5°  for  N = 4 , c = .8. 
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List  of  Figures 


Each  of  the  figures  A.  78  - A.  83  consists  of  5 sub- 
figures labeled  a,b,c,d,e.  The  subfigures  are  plots 
of  the  following  quantities: 

a:  e^t) 
b:  02 ( t) 
c:  63 (t) 
d:  04(t) 

e:  Dk(t)  , k = 1,2, 3, 4. 

For  example,  Figure  A.  72  c represents  the  component 
03vt)  . 


Figure 

Number 

Colatitude 
of  Vortices 

Coriolis  Parameter 

Initial 

Displacement 

A. 78 

35° 

.8 

= 0.5° 

A.  79 

35° 

.8 

e 8X  = -15.0° 

A. 80 

35° 

.8 

za1  = -15  0 

A. 81 

42° 

.8 

1 

c$1  = .25° 

A. 82 

42° 

.8 

ea^  = .25° 

A. 83 

50° 

.8 

e&1  = 0.1  0 

2 
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Figure  A. 78  d 
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Figure  A. 78  e 
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